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AN  ASSESSMENT  OF 
RESERVOIR  DENSITY  CURRENTS 
AND  INFLOW  PROCESSES 

PART  I:  INTRODUCTION 

1.  When  a  reservoir  is  formed  by  impounding  a  free-flowing  river, 
both  the  flow  regime  and  ecosystem  characteristics  change.  The  type  of 
ecosystem  that  develops,  the  resultant  water  quality,  and  the  ability  to 
meet  project  purposes  is  dependent  on  understanding  how  inflows  enter 
and  pass  through  impoundments.  Since  tributary  inflows  are  usually  the 
major  sources  of  materials  such  as  nutrient  and  suspended  and  dissolved 
matter  which  impact  reservoir  water  quality,  techniques  for  predicting 
the  characteristics  of  density  currents  are  required  to  determine  the 
ways  in  which  management  and  operational  alternatives  can  regulate 
reservoir  and  release  water  quality. 

2.  The  objectives  of  this  study  were  to: 

■a.  Review  the  literature  and  field  studies  on  density  currents 
in  reservoirs. 

b^.  Develop  techniques  for  predicting  the  placement,  plunge 
point  location,  thickness,  propagation  speed,  and  mixing 
of  density  currents. 

3.  The  predictive  techniques  include  both  simple  computational 
equations  (e.g.,  to  calculate  how  fast  a  density  current  will  move 
through  a  reservoir  under  the  influence  of  a  specific  set  of  circum¬ 
stances)  and  formalized  algorithms  for  complex  water  quality  models 
(e.g.,  CE-QUAL-R1  (Environmental  Laboratory  1982)).  The  predictive 
techniques  must,  therefore,  be  sufficiently  general  to  encompass  the 
many  different  reservoir  types  operated  by  the  U.  S.  Army  Corps  of 
Engineers  (CE)  and  also  be  formulated  on  data  and  measurements  routinely 
available  for  CE  reservoirs. 
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PART  II:  BACKGROUND 


4.  Since  the  tributary  inflow  density  usually  differs  from  the 
density  of  the  reservoir  surface  water,  inflows  enter  and  move  through 
reservoirs  as  density  currents.  Bell  (1942)  defines  a  density  current 
as  "a  gravity  flow  of  a  liquid  or  a  gas  through,  under,  or  over  a  fluid 
of  approximately  equal  density."  Density  flows  differ  from  normal  fluid 
flows  because  the  buoyancy  of  the  surrounding  fluid  reduces  the  gravity 
force  by  the  normalized  density  difference,  Ap/p  (i.e.  ,  reduced 
gravity  force  =  ^g  ).  Vertical  movements  such  as  waves  are  therefore 
greater  at  the  density  interface  than  at  the  air-water  interface. 
Katabatic  winds,  dust  storms  (e.g.,  Sudan  haboob) ,  sea  breezes,  and  the 
advancement  of  a  cold  front  are  other  examples  of  naturally  occurring 
density  currents. 

5.  In  reservoirs,  density  currents  are  caused  by  density  dif¬ 
ferences  due  to  temperature,  total  dissolved  solids  (TDS) ,  and  suspended 
solids  (SS).  In  most  reservoirs,  however,  density  differences  are 
predominantly  caused  by  temperature  differences.  For  example,  at  25°C 
it  takes  approximately  330  mg/1  of  dissolved  solids  or  420  mg/1  of 
suspended  solids  (specific  gravity  *  2.65)  to  equal  the  density  differ¬ 
ence  caused  by  a  1°C  temperature  change.  Procedures  to  compute  water 
densities  based  on  temperature,  TDS,  and  SS  concentrations  are  presented 
in  Appendix  A. 


Qualitative  Description 

6.  When  a  tributary  inflow  enters  a  reservoir,  it  displaces  the 
reservoir  water  ahead  of  it.  Increases  in  water  depth  and  cross- 
sectional  area  cause  flow  velocities  to  decrease.  If  there  is  no 
density  difference  between  the  tributary  inflow  and  reservoir  waters, 
the  point  of  maximum  velocity  remains  near  the  water  surface  as  in  open 
channel  flow.  The  major  forces  acting  on  a  parcel  of  water  include  the 
dynamic  force  exerted  by  the  stream  momentum  and  resisting  force  due  to 
bed  shear.  Wind  shear  at  the  air-water  interface  can  be  either  a 


positive  or  negative  force.  As  the  inflowing  parcel  of  water  moves 
towards  the  dam,  concentrations  of  inflowing  constituents  decrease  due 
to  diffusion  and  dispersion. 

7.  In  most  instances,  the  tributary  inflows  differ  in  density 
from  the  reservoir  waters.  Depending  on  the  magnitude  of  this  den¬ 
sity  difference,  density  currents  can  enter  the  epilimnion,  metalim- 
nion,  or  hypolimnion  (Figure  1).  When  the  inflow  density  is  less  than 
the  surface  water  density,  the  inflow  will  float  on  the  water  surface 
(i.e.,  an  overflow).  In  an  overflow,  the  excess  hydrostatic  pressure  in 
the  density  current  causes  the  current  to  flow  in  all  directions  not 
obstructed  by  boundaries.  Overflows  are  susceptible  to  mixing  from  wind- 
induced  mechanisms  and  diurnal  heating  and  cooling  processes.  Wind 
shear  can  direct  the  overflow  into  a  cove  or  prevent  it  from  moving 
downstream.  Vertical  mixing  can  also  weaken  the  density  difference  by 
mixing  the  overflow  throughout  the  water  column,  causing  it  to  lose  its 
integrity. 

8.  If  the  inflow  density  is  greater  than  the  water  surface  density, 
the  inflow  will  push  the  reservoir  water  ahead  until  the  buoyancy  forces 
dominate  and  the  inflow  plunges  beneath  the  water  surface.  The  plunge 
point  is  sometimes  made  visible  because  of  turbidity  or  the  accumulation 
of  floating  debris,  which  may  indicate  a  stagnation  point.  The  location 
of  the  plunge  point  is  determined  by  a  balance  between  the  stream  mo¬ 
mentum,  the  pressure  gradient  across  the  interface  separating  the  river 
and  reservoir  waters,  and  the  resisting  shear  forces.  Some  mixing 
(termed  initial  mixing)  occurs  at  the  plunge  point  because  of  the  large 
eddies  formed  by  flow  reversals  and  pooling  of  the  inflowing  water 
(Akiyama  and  Stefan  1981).  Knapp  (1942)  noticed  that  the  flow  in  the 
vicinity  of  the  plunge  point  occurred  at  the  bottom  of  this  pooled  mixing 
zone  (Figure  2).  Ford  et  al.  (1980)  and  Kennedy  et  al.  (1983)  substanti¬ 
ated  this  pooling  phenomenon  during  dye  studies  on  DeGray  Lake,  Arkansas, 
and  West  Point  Reservoir,  Georgia,  when  the  dye  clouds  appeared  to  have 
stalled  at  the  plunge  point  (Figure  3).  The  location  of  the  plunge 
point  can  also  be  influenced  by  morphological  factors.  Changes  in  the 
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bed  slope  (e.g.,  due  to  sediment  deposition),  bed  friction,  and  cross- 
sectional  area  may  affect  its  location.  For  a  river  entering  a  wide 
lake,  the  plunge  point  may  actually  be  a  point  rather  than  a  line. 

9.  After  the  inflow  plunges,  it  follows  the  old  river  channel 
(thalweg)  as  an  underflow.  The  speed  and  thickness  of  the  underflow 

is  determined  by  a  flow  balance  between  the  shear  forces  and  the  accel¬ 
eration  due  to  gravity  (i.e.,  gradually  varying  flow  theory).  The 
momentum  equation  used  in  two-layered  flow  analysis  is  similar  to 
the  open-channel  flow  formulation  except  for  the  additional  buoyancy 
force.  An  underflow  will  entrain  overlying  reservoir  water  due  to 
turbulence  generated  by  bottom  roughness.  Changes  in  the  underflow  den¬ 
sity  from  entrainment  must  be  quantified  before  a  density  interflow  or 
intrusion  can  be  analyzed,  or  the  vertical  placement  of  an  interflow 
will  be  incorrect. 

10.  An  interflow  or  intrusion  occurs  when  a  density  current  leaves 
the  river  bottom  and  propagates  horizontally  into  a  stratified  body  of 
water.  Intrusions  differ  from  overflows  and  underflows  because  an  in¬ 
trusion  moves  through  a  reservoir  at  a  elevation  where  the  intrusion 
and  reservoir  densities  are  similar.  Intrusions  require  a  continuous 
inflow  and/or  outflow  for  movement,  or  they  stall  and  collapse  (i.e., 
dissipate) .  Entrainment  is  usually  neglected  in  the  analysis  of  an 
intrusion  since  the  density  gradient  in  the  metalimnion  creates  strong 
buoyancy  forces  which  inhibit  mixing. 

Field  Observations 

11.  The  existence  of  density  currents  in  reservoirs  has  been  known 
for  many  years.  Early  studies  consisted  of  watching  turbid  waters  flow 
into  a  reservoir  and  disappear.  Observations  of  turbid  release  water 
from  Elephant  Butte  Dam,  Rio  Grande  River,  New  Mexico,  in  the  early 
1900's  indicated  that  turbid  inflows  could  flow  through  a  reservoir  in  a 
distinct  layer  (Grover  and  Howard  1938) .  Grover  and  Howard  (1938)  and 
accompanying  discussions  described  the  passage  of  silt-laden  Colorado 
River  water  through  Lake  Mead,  Nevada,  based  on  chemical  analysis  of 
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inflow  and  outflow  water.  They  suggested  the  possibility  of  increasing 
the  life  of  a  reservoir  by  promoting  the  passage  of  fine  silt  through 
reservoirs  by  bottom  withdrawal.  Uiebe  (1939)  suggested  the  existence 
of  density  currents  in  Norris  Reservoir,  Tennessee,  based  on  routine 
limnological  measurements  including  dissolved  oxygen,  alkalinity,  pH, 
nitrite  nitrogen,  and  turbidity. 

12.  Recent  observations  of  density  currents  have  been  more  quanti¬ 
tative  using  both  naturally  occurring  substances  and  fluorescent  dyes  to 
trace  density  currents.  Direct  measurements  of  density  currents  have 
been  limited  since  the  velocities  are  below  the  detection  limits  of  most 
current  meters.  The  exception  was  the  velocity  measurements  made  by  the 
Tennessee  Valley  Authority  (TVA)  using  the  DWICA  isotopic  current 
analyzer.  These  measurements  were,  however,  time-consuming  and  there¬ 
fore  limited  (Elder  and  Wunderlich  1972) .  Unsteady  flow  also  made  data 
analyses  difficult. 

Natural  tracers 

13.  In  many  reservoirs,  the  properties  of  the  inflowing  streams 
are  sufficiently  different  from  the  reservoir  waters  that  the  inflows 
can  be  tracked  through  the  reservoir.  Differences  in  turbidity,  TDS, 
specific  conductance,  temperature,  and  other  naturally  occurring  sub¬ 
stances  have  all  been  used  to  trace  inflows.  Fry  et  al.  (1953)  docu¬ 
mented  density  currents  in  17  TVA  reservoirs  based  on  temperature, 
alkalinity,  and  chloride  concentrations.  The  impact  of  these  currents 
on  industrial  water  use,  stream  pollution,  recreation,  municipal  water 
supply,  and  steam  plant  condensing  water  was  described  in  detail.  In 
many  cases,  the  density  currents  caused  reverse  currents  which  moved 
undesirable  contaminants  upstream  into  water  intakes.  Other  studies  of 
density  currents  in  TVA  reservoirs  (e.g.,  Churchill  1958,  Churchill 
1965,  Elder  and  Wunderlicl.  1972)  described  similar  experiences.  Other 
observations  of  density  currents  include  Bighorn  Lake,  Montana  (Soltero 
et  al.  1974),  Lake  Powell,  Utah-Arizona  (Johnson  and  Merritt  1979),  and 
Lake  Kootenay,  British  Columbia  (Hamblin  and  Carmack  1978,  Carmack  1979), 
among  others. 
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14.  Serruya  (1974)  described  the  inflow  into  Lake  Kinneret,  Israel, 
using  conductivity  data.  This  was  an  excellent  example  of  inflow  into  a 
wide  body  of  water.  It  also  illustrated  the  effects  of  Coriolis  accel¬ 
erations  which  caused  the  density  interfaces  to  tilt  to  the  right  side 
(looking  downstream)  of  the  basin,  as  expected  in  the  northern  hemis¬ 
phere.  Hamblin  and  Carmack  (1978)  have  analyzed  these  types  of  flow  in 
detail. 

15.  Hebbert  et  al.  (1979)  studied  a  saline  underflow  into  Welling¬ 
ton  Reservoir,  Australia,  where  the  flushing  of  salts  from  the  watershed 
caused  salinity  concentrations  to  peak  before  flow.  The  authors  showed 
that  both  the  flood  and  salinity  surges  travelled  at  constant  speeds  but 
that  the  flood  surge  travelled  2.3  times  faster  than  the  salinity  surge 
and  therefore  moved  past  the  salinity  surge  in  the  reservoir. 

16.  Perhaps  the  most  complete  set  of  data  describing  turbid 
inflows  into  a  reservoir  has  been  collected  by  J.  Nix  and  staff, 

Ouachita  Baptist  University,  on  UeGray  Lake,  Arkansas.  Observations 
included  overflows,  interflows,  and  underflows  during  both  base  and 
elevated  flow  regimes.  Inflows  were  tracked  using  turbidity,  calcium, 
and  specific  conductance.  These  density  currents  are  discussed  in  part 
by  Thornton  et  al.  (1980),  Ford  et  al.  (1980),  Ford  and  Johnson  (1981), 
and  Nix  (1981)  and  completely  documented  by  Johnson  et  al.  (1983). 

Dye  studies 

17.  When  there  is  no  naturally  occurring  property  which  can  be 
easily  detected  and  measured,  inflows  can  be  labeled  with  an  artificial 
tracer  to  follow  their  behavior.  Fluorescent  dyes  have  proven  to  be 
excellent  tracers  because  of  their  properties  of  high  detectibility  and 
ease  of  handling.  Johnson  (1983)  discusses  fluorescent  dyes  and  their 
use. 

18.  A  number  of  field  studies  using  fluorescent  dyes  have  been 
performed.  Many  of  these  studies  have  not  been  fully  documented  or 
published  in  a  form  that  is  readily  available,  thus,  the  following 
discussion  is  not  complete. 
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19.  One  of  the  first  dye  studies  documented  was  by  Elder  and 
Wunderlich  (1968) .  They  dyed  the  Nantahala  River  flowing  into  Fontana 
Reservoir,  Tennessee,  using  Rhodamine  B  fluorescent  dye.  The  dye  flowed 
through  the  reservoir  as  an  interflow  and  travelled  approximately  32  km 
in  a  period  of  6  weeks  (approximately  0.9  cm/s)  before  becoming  unde¬ 
tect  ible.  As  the  reservoir  elevation  dropped,  the  dye  also  sank  to  a 
lower  elevation  but  remained  well-defined  indicating  that  little  or  no 
mixing  occurred.  Since  the  inflow  temperature  varied  by  over  2°C 
during  dye  injection,  it  was  suspected  that  the  dye  moved  into  dif¬ 
ferent  layers,  creating  the  appearance  of  vertical  mixing  or  diffusion. 

20.  Bayne  (1967)  studied  the  movement  and  dispersion  of  the 
Chattahoochee  River  as  it  moved  through  Lake  Eufaula,  Alabama,  using 
Rhodamine  WT  fluorescent  dye.  In  the  upper  reaches  of  the  reservoir, 
the  water  movement  was  riverlike  with  an  average  velocity  of  0.13  m/s. 

In  the  middle  section,  the  dye  moved  as  an  underflow  following  the  old 
river  channel  with  an  average  speed  of  0.04  m/s.  In  the  lower  end  of 
the  reservoir,  the  dye  movement  was  erratic  and  was  even  observed  to 
move  upstream.  This  may  have  been  caused  by  a  constriction  in  an 
earthen  causeway. 

21.  Morris  and  Thackston  (1969)  used  two  dye  injections  to  study 
the  movement  of  the  South  and  North  Forks  of  the  Holston  River  through 
Cherokee  Reservoir,  Tennessee.  In  the  first  dye  study,  a  cold-water 
pulse  from  the  South  Fork  was  dyed;  in  the  second  study,  the  warmer 
North  Fork  was  dyed.  Both  inflows  moved  through  the  reservoir  as  inter¬ 
flows;  during  the  second  study,  the  interflow  was  thicker  and  proceeded 
at  a  greater  depth  because  stratification  was  weaker.  In  both  studies, 
the  dye  was  well-mixed  laterally.  During  a  period  of  slow  continuous 
drawdown,  no  dye  was  observed  to  move  into  the  coves.  Some  horizontal 
movement  of  dye  into  the  coves  was  observed  when  the  water  surface 
fluctuated  by  up  to  0.12  m.  Virtually  no  vertical  transport  was 
observed. 

22.  Chasse  and  Slotta  (1972)  made  12  dye  releases  into  three 
reservoirs  (Hills  Creek,  Cougar,  and  Dorena)  in  Oregon.  Although  most 


of  the  studies  were  on  Hills  Creek  Reservoir,  a  wide  range  of  conditions 
was  studied.  Both  interflows  and  underflows  were  observed;  current 
thicknesses  varied  from  3  to  12  m;  current  speeds  varied  from  0.02  to 
0.13  m/s.  Chasse  and  Slotta  also  computed  entrainment  rates,  which 
varied  from  1.00  (no  entrainment)  to  4.0  (300%  entrainment).  In  Cougar 
Reservoir,  the  dye  entered  at  two  elevations;  no  explanation  was  given 
for  the  formation  of  two  distinct  interflows. 

23.  Sutron  Corporation  (1980)  studied  the  lateral  mixing  charac¬ 
teristics  in  Lewis  and  Clark  Lake,  South  Dakota,  on  the  Missouri  River. 
Rhodamine  WT  fluorescent  dye  was  injected  into  the  river  above  the  lake 
and  followed  through  the  lake.  In  the  upper  reaches  of  the  lake,  the 
dye  remained  in  the  zone  of  maximum  velocity.  In  the  main  body  of  the 
lake,  wind  appeared  to  influence  the  movement  of  the  dye.  The  dye  was 
not  well-mixed  laterally  at  any  time.  A  two-dimensional  transport  model 
was  used  to  investigate  lateral  mixing  and  the  influence  of  the  wind. 

24.  Fischer  and  Smith  (1982)  studied  an  underflow  in  the  Las  Vegas 
arm  of  Lake  Mead.  Rhodamine  WT  dye  was  used  to  determine  the  movement 
of  inflowing  nutrients  into  the  water  surface.  The  study  demonstrated 
that  the  nutrients  (dye)  reached  the  surface  through  the  current  field 
generated  by  large  internal  waves  (as  high  as  6  m)  and  the  transient 
formation  and  destruction  of  secondary  thermoclines  due  to  variations  in 
meteorological  conditions. 

25.  Ford  et  al.  (1980)  also  documented  the  effects  of  variations 
in  meteorological  conditions  on  inflow  currents.  Convective  mixing  at 
night  divided  a  dye  cloud  into  three  distinct  parts.  This  study  and 
others  on  DeGray  Lake,  described  by  Johnson  et  al.  (1983),  showed  that 
even  under  low-flow  conditions,  density  currents  move  through  lakes  in 
distinct  layers. 

Complicating  Factors 

26.  Based  on  the  field  observations,  there  are  a  number  of  factors 
which  complicate  the  analysis  and  prediction  of  density  currents  in 
reservoirs.  These  include  complicated  morphometry,  unsteady  and 
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Figure  4.  Typical  long,  dendritic  reservoir  with  meandering  thalweg. 

multiple  tributary  inflows,  variable  stratification,  and  man-made 
perturbations  and  constraints. 

Complicated  morphometry 

27.  Reservoirs  are  typically  built  in  deep  valleys  downstream  of 
the  confluence  of  several  tributaries.  The  resulting  lake  is  long, 
highly  dendritic,  and  possibly  surrounded  by  high  terrain  (e.g.  Figure 
4).  The  many  coves,  embayments,  and  islands  make  the  width  highly 
variable  and  difficult  to  define.  Field  observations  of  turbidity 
currents  in  DeGray  Lake,  Arkansas,  (Ford  and  Johnson  1981)  indicated 
that  in  many  instances,  density  currents  were  not  well-mixed  laterally 
but  followed  the  thalweg  of  the  inundated  river  (Figure  5).  This  phe¬ 
nomenon  has  also  been  observed  in  other  reservoirs  (e.g.  Lake  Red  Rock, 
Iowa,  Lewis  and  Clark  Lake,  South  Dakota).  In  many  cases,  especially 
in  reservoirs  that  have  not  been  cleared  of  trees,  the  zone  above  the 
old  channel  may  be  the  path  of  least  resistance.  Sonar  recordings  from 
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Downstream  of  Station  10,  the  trees  were  topped  at  elevation  113  m  above 
mean  sea  level  (msl) ,  or  approximately  11  m  below  normal  pool  elevation 
(Figure  6b),  an  elevation  within  the  metalimnion  during  the  summer 
stratification.  The  cross  section  shown  in  Figure  6c  is  located  near 
Station  4.  The  old  river  channel  is  well-defined  in  all  three  cross 
sections.  Density  flows  occurring  at  elevations  below  the  tree  level 
were  constrained  to  the  old  river  channel;  many  of  the  interflows 
occurred  at  elevations  near  or  slightly  above  the  tree  tops. 

28.  Since  density  currents  may  not  be  well-mixed  laterally,  it 
may  be  necessary  to  define  the  width  or  limits  for  the  zone  of  convey¬ 
ance.  The  appropriate  width  and  length  scales  for  density  currents  may 
be  the  width  and  length  of  the  old  river  channel,  respectively. 

29.  Old  roadways,  railroad  crossings,  etc.,  that  were  not  removed 
prior  to  the  filling  of  a  reservoir  may  also  act  as  hydraulic  controls 
to  direct  the  density  current.  For  example,  in  Lake  Red  Rock  the 
embankment  for  a  low-water  bridge  was  not  removed  (Figure  7).  This 
embankment,  which  is  submerged  at  flood  pool,  acts  as  a  barrier  to 
direct  flow  along  the  south  shore.  Old  cofferdams,  construction  roads, 
etc.,  are  also  sometimes  inundated  causing  perturbations  and  irregu¬ 
larities  in  the  anticipated  flow  field.  If  an  old  roadway  extends 
across  the  reservoir,  the  density  current  may  buildup  behind  the  road¬ 
way  prior  to  flowing  over  it  (Knapp  1942). 

30.  The  slope  of  old  river  channel  may  be  modified  due  to  sedimen¬ 
tation  and  changes  in  the  configuration  of  the  delta  region.  In  Lake 
Red  Rock,  over  6  m  of  sedimentation  occurred  in  the  old  river  channel 
during  the  first  seven  years  of  operation  (Kennedy  et  al.  1981).  The 
path  of  the  channel  through  the  delta  also  changed  several  times.  In 
Lake  Mead,  Nevada-Arizona,  the  delta  progressed  67  km  in  13  years 
(Vanoni  1975). 

Unsteady  and  multiple  inflows 

31.  The  magnitude  and  timing  of  tributary  inflows  to  reservoirs  is 
dependent  on  the  size  and  shape  of  the  watershed,  antecedent  conditions, 
and  the  distribution  of  precipitation.  In  northern  latitudes,  runoff 
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from  the  spring  snowmelt  may  dominate  the  water  budget;  elsewhere,  the 
water  budget  may  be  dominated  by  individual  storm  events  (e.g.,  see 
Figure  8).  In  either  case,  the  major  contributions  to  a  reservoir's 
water  budget  occur  during  dynamic  elevated  flow  events  and  not  during 
steady  base  flow  conditions. 

32.  The  density  of  the  inflowing  waters  is  also  variable.  As 
previously  stated,  water  density  is  a  function  of  temperature,  sus¬ 
pended  solids  (SS),  and  total  dissolved  solids,  (TDS) .  Water  tempera¬ 
tures  vary  seasonally,  synoptically  (i.e.,  for  periods  of  5-10  days), 
and  diurnally  in  response  to  hydrometeorological  forcing  (Ford  and  Ford 
1982).  SS  and  TDS  concentrations  vary  with  runoff,  flow,  antecedent 
conditions,  and  other  factors.  Examples  of  the  types  of  variations  in 
flow,  temperature,  SS,  TDS,  and  density  that  can  occur  during  a  storm 
event  are  shown  in  Figures  9  through  11  for  DeGray  Lake.  Flows  in- 

3 

creased  from  approximately  8.5  m  /s  at  1200  hours  on  3  May  to  a  peak  of 

3 

190  m  /s  at  0700  hours  on  4  May,  while  stream  temperatures  dropped  from 
18. 2°C  on  3  May  to  13.8°C  on  5  May.  Turbidity  increased  on  the  rising 
limb  of  the  hydrograph  and  attained  a  maximum  of  73  NTU's  11  hours  prior 
to  the  peak  of  the  hydrograph.  In  contrast,  specific  conductance,  which 
is  proportional  to  TDS,  decreased  with  increased  flow.  Since  water 
density  is  dependent  on  temperature,  TDS,  and  SS,  it  will  also  vary 
significantly  during  the  storm  event.  As  the  inflow  density  varies 
during  the  storm,  the  inflow  will  be  placed  at  different  levels  in  the 
reservoir  (Figure  12).  Other  storm  events,  described  by  Thornton  et  al. 
(1980),  Ford  and  Johnson  (1981),  and  Johnson  et  al.  (1983),  demonstrated 
that  each  storm  is  different  and,  in  some  cases,  temperature  increases 
during  the  storm  (Figure  13).  Runoff  also  can  increase  TDS  or  salinity 
on  the  rising  side  of  the  hydrograph  (e.g.,  Hebbert  et  al.  1979). 

33.  The  unique  thermal-density  properties  of  water  also  create 
some  complicated  inflow  patterns  as  water  warms  and  passes  through  its 
maximum  density  at  4°C.  Using  data  from  Kamloops  Lake,  British  Columbia, 
Carmack  (1979)  found  that  when  inflow  water  above  4°C  enters  a  lake  with 
a  temperature  of  less  than  4°C,  the  two  waters  mix  creating  a  denser 
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Figure  11.  Temporal  variation  of  inflow  and  inflow  density  during  June  1977 
storm  event,  Caddo  River,  DeGray  Lake,  AR. 
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Figure  14.  Inflow  mixing  and  cabbeling  as  inflow  temperature 
increases  from  0°C  to  8°C  (after  Carmack  1979). 

mixture  near  4°C,  which  sinks.  This  process  is  called  cabbeling. 
Following  Carmack  (1979),  a  simple  example  which  assumes  a  constant  lake 
temperature  of  2°C  and  an  inflow  temperature  that  warms  from  0  to  »°C, 
is  used  to  illustrate  the  complex  circulation  patterns  which  can  develop 
(Figure  14).  In  Figure  14a,  the  inflowing  water  at  a  temperature  of 
0°C  is  less  dense  than  the  lake  water  at  2°C  and  the  inflow  enters  as  su 
overflow.  In  Figure  14b,  the  inflow  is  more  dense  and  enters  the  lake 
as  an  underflow.  When  the  river  temperature  is  6°C,  it  mixes  with  the 
2°C  lake  water  to  form  a  mixture  which  is  more  dense  and  sinks  (Figure 
14c).  If  the  river  temperature  is  8°C,  it  may  form  a  mixture  with  the 
lake  water  (a)  in  the  range  of  2  to  6°C  which  would  sink  and  create 
mixing  or  (b)  in  the  range  of  6  to  8°C,  in  which  case  the  inflow  would 
continue  to  float  on  the  water  surface. 

34.  Instabilities  in  the  inflow  current  can  also  occur  when  water 
density  is  dependent  on  SS  and  TDS  concentrations.  As  velocities  de¬ 
crease  when  the  current  enters  a  reservoir,  SS  will  settle  out  thereby 
changing  the  water  density.  Since  dissolved  salts  and  temperature 
diffuse  at  different  rates  (molecular  diffusivity  of  heat  is  greater 


chan  Chat  of  salt)  instabilities  can  also  occur  which  increase  mixing. 
This  process  is  called  double  diffusion  and  is  described  by  Turner 
(1973). 

35.  When  multiple  tributaries  enter  a  reservoir,  they  may  move 
through  as  distinct  currents  or  they  may  interact  to  form  a  common 
current.  Koryak  et  al.  (1979)  found  two  distinct  currents  in  Tygart 
River  Lake,  West  Virginia.  The  smaller,  cooler  tributary  flowed  through 
the  lake  as  a  distinct  interflow  at  a  deeper  stratum.  In  West  Point 
Reservoir,  Georgia,  the  warmer,  clearer  waters  of  Yellow  Jacket  Creek 
appear  to  flow  over  the  heavier,  sediment-laden  waters  of  the  Chat¬ 
tahoochee  River  (Kennedy  et  al.  1982).  The  interactions  between  these 
two  tributaries  are  highly  complex. 

Variable  Stratification 

36.  As  with  the  tributary  inflows,  the  in-lake  stratification  also 
varies  in  response  to  hydrometeorologic  forcing,  but  not  at  the  same 
rate.  The  larger  volume  of  water  in  the  lake  responds  at  a  slower  rate 
than  the  river  water.  During  the  spring,  the  river  water  warms  at  a 
faster  rate  than  the  lake  water;  in  the  fall,  river  water  cools  at  a 
faster  rate  than  the  lake  water.  This  is  why  overflows  typically  occur 
in  the  spring  and  interflows  and  underflows  occur  in  late  summer  and 
fall  (Figure  8) . 

37.  An  example  of  changes  in  stratification  occurring  during  a 
storm  event  at  DeGray  Lake  is  shown  in  Figure  15.  An  intrusion  entered 
the  reservoir  at  an  elevation  of  116.5  m  msl,  or  at  a  temperature  of 
20°C  as  indicated  by  the  intersection  of  the  temperature  profiles  in 
Figure  15.  On  18  June,  the  temperature  stratification  was  altered  due 
to  the  intrusion.  The  turbidity  and  temperature  profiles  for  that  day 
indicated  that  the  intrusion  cooled  the  water  above  and  warmed  the  water 
below  the  entrance  elevation.  By  19  June,  the  water  that  was  warmed 
below  the  center-line  elevation  had  returned  to  its  previous  temperature, 
indicating  that  the  stratification  was  only  temporarily  displaced  on  18 
June  and  that  no  significant  mixing  had  occurred.  In  contrast,  the 
stratification  above  the  center-line  elevation  was  permanently  altered. 
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Figure  15.  Variations  in  thermal  stratification  during 
June  1977  storm  event,  DeGray  Lake,  AR. 

indicating  mixing.  This  mixing  was  probably  due  to  a  shear  mechanism. 
After  19  June,  the  temperature  profiles  showed  that  surface  mixing 
dominated. 

38.  Surface  mixing  phenomena  were  observed  to  dominate  during  a 
dye  study  on  DeGray  Lake  in  1979.  Nocturnal  convective  mixing  caused 
the  dye  cloud  to  break  into  two  parts  on  14  October  and  into  three  parts 
on  15  October  (Figure  16).  The  origin  of  these  breaks  can  be  traced 
back  to  the  locations  at  kilometers  2.8  and  4.0  where  a  greater  volume  of 
water  from  deep  holes  and  embayments  was  available  for  mixing  and  dilu¬ 
tion  of  the  dye  cloud  at  night.  This  convective  mixing  contributed  to 
the  increased  diffusion  coefficient.  It  was  only  after  the  weather  had 
warmed  and  the  reservoir  had  restratified  that  the  classical  underflow 
was  observed  on  17  October. 
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Figure  16.  Nocturnal  mixing  during  October 


Man-made  perturbations  and  constraints 


39.  Project  operation  and  point  source  discharges  also  can  modify 
the  behavior  of  density  currents.  Many  reservoirs  are  characterized  by 
wide  fluctuations  in  pool  elevation.  Water  is  stored  during  periods  of 
high  flow  and  released  during  periods  of  low  flow.  In  hydropower  pro¬ 
jects,  releases  are  made  only  to  generate  power.  In  either  case,  when¬ 
ever  a  project  is  operated  such  that  inflow  is  not  equal  to  outflow, 
density  currents  are  governed  by  unsteady  forces  making  their  analyses 
difficult.  For  example.  Ford  and  Johnson  (1981)  showed  that  a  density 
flow  stalled  and  collapsed  in  DeGray  Lake  when  the  outflows  were  stopped. 
In  both  the  May  1979  storm  (Figure  9)  and  October  1976  storm  (Figure 
13),  the  lake  surface  rose  since  the  inflow  was  greater  than  the  outflow. 
This  caused  water  to  move  back  into  the  coves,  complicating  the  analysis 
of  the  density  current. 

40.  The  level  of  the  withdrawal  is  also  important.  Bottom  with¬ 
drawal  projects  reinforce  the  movement  of  underflows.  Withdrawing  water 
from  a  level  different  from  that  of  the  density  current  can  generate 
secondary  currents  which  complicate  the  analysis. 

Water  Quality  Implications 

41.  Inflows  to  impoundments  not  only  contribute  to  the  water 
budget  but  also  carry  solids,  nutrients,  bacteria,  and  other  nubstances 
that  affect  reservoir  water  quality.  These  materials  may  enter  the 
epilimnion,  metalimnion,  or  hypoliranion  depending  on  whether  an  overflow, 
interflow,  or  underflow  occurs.  If  an  overflow  occurs,  the  nutrients 
enter  the  epilimnion  directly  and  are  immediately  available  for  phyto¬ 
plankton  uptake.  Coliform  bacteria  entering  the  epilimnion  may  impact 
body-contact  recreation  in  the  upper  reaches  of  a  reservoir  (Thornton, 

et  al.  1980).  If  an  interflow  or  underflow  occurs,  the  transported 
materials  will  plunge  below  the  water  surface  and  may  pass  through  the 
reservoir  at  a  level  that  would  not  impact  the  surface  waters.  Fischer 
and  Smith  (1982)  showed  that  it  is  still  possible  for  inflowing  nutri¬ 
ents  to  get  to  the  surface  waters  by  various  mixing  mechanisms  even 
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though  a  well  defined  interflow  existed.  A  cold  inflow,  rich  in  dis¬ 
solved  oxygen,  may  also  find  its  way  into  an  anaerobic  hypolimnion, 
allowing  oxidation  of  substances  there.  Koryak  et  al.  (1979)  showed 
that  interflows  in  Tygart  River  Lake,  West  Virginia,  helped  to  mitigate 
the  effects  of  acid  mine  drainage. 

42.  Recent  studies  have  shown  that  significantly  more  material 
enters  a  reservoir  during  storm  events  than  during  base  flow  conditions 
(Perrier  et  al.  1977,  Westerdahl  et  al.  1981).  During  storms,  certain 
constituents  (e.g.,  phosphorus,  coliform  bacteria,  turbidity''  charac¬ 
teristically  load  on  the  rising  side  of  the  hydrograph  while  others 
(e.g.,  nitrate  and  many  metals)  load  on  the  falling  limb  (Perrier  et  al . 
1977,  Westerdahl  et  al.  1981).  The  ultimate  fate  of  a  constituent  will 
depend  not  only  on  how  it  moves  through  the  reservoir  but  also  on  how  it 
loads  on  the  hydrograph.  Using  the  two-dimensional  model  LARM,  Johnson 
et  al.  (1981)  showed  that  turbidity,  fecal  coliforms,  nitrogen,  and 
phosphorus  moved  through  DeGray  Lake  at  different  rates  and  elevations 
depending  on  the  unique  manner  in  which  each  constituent  loaded  on  the 
hydrograph. 

43.  Density  currents  can  also  produce  flow  reversals  in  the  reser¬ 
voir  that  cause  undesirable  materials  to  be  transported  upstream.  When 
the  cool  inflow  plunges  beneath  the  surface  and  enters  the  reservoir  as 
an  interflow  or  as  an  underflow,  some  of  the  reservoir  water  is  en¬ 
trained  and  a  counteracting  circulation  pattern  is  set  up  in  the  over- 
lying  surface  waters.  This  reverse  flow  pattern  transports  floating 
debris  upstream  to  the  plunge  point  where  it  remains  at  the  surface. 
Water  is  again  entrained  and  plunges  to  continue  the  circular  pattern. 
Bell  (1942)  noted  that  in  June  1941  a  driftwood  barrier  was  set  up  in 
Lake  Mead,  Nevada-Arizona,  in  this  manner  from  bank  to  bank,  extending 
several  thousand  feet  long  and  effectively  blocking  all  navigation. 

44.  Another  problem  associated  with  induced  flow  reversal  occurs 
when'xmunicipal  water  intakes  are  located  upstream  of  waste  discharges. 
Fry  et  al.  (1953)  reported  two  such  cases  for  the  cities  of  Harriman  and 
Knoxville,  Tennessee.  Knoxville's  water  intakes  are  located  just 


inside  the  upstream  limit  of  Fort  Loudin  Reservoir's  pool,  upstream  of 
sewage  and  industrial  waste  outfalls.  A  TVA  study  determined  that, 
unless  the  tributary  inflows  to  the  reservoir  filled  the  entire  river¬ 
like  channel  in  the  reservoir  near  the  intakes,  the  inflow  would  plunge 
and  become  an  underflow.  Then,  as  the  underflow  moved  downstream,  the 
resulting  return  current  above  it  would  transport  sewage  and  industrial 
wastes  upstream  to  Knoxville's  water  intakes.  To  prevent  such  an 
occurrence,  minimum  release  criteria  were  established  for  the  two  up¬ 
stream  reservoirs  supplying  the  cold  tributary  inflows  to  Fort  Loudin 
Reservoir. 

45.  At  Harriman  a  more  complicated  situation  existed.  Harriman's 
water  intakes  were  located  in  the  Emory  River  arm  of  Watts  Bar  Reser¬ 
voir  upstream  of  sewage  and  paper  mill  discharges.  Density  currents 

in  the  Clinch  River  arm  of  the  reservoir  often  traveled  downstream  into 
a  pool  connecting  both  arms  and  then  upstream  into  the  Emory  River  arm. 
The  invasion  of  a  density  current  into  the  Emory  River  arm  itself  posed 
no  hazard  to  Harriman.  A  problem  was  created,  however,  when  the  in¬ 
truding  density  current  stalled  and  began  moving  back  downstream  the 
Emory  River  arm  as  a  result  of  reductions  in  outflows  from  the  upstream 
reservoir  feeding  the  density  current.  As  the  density  current  moved 
back  downstream,  downstream  water  moved  upstream  to  replace  it,  trans¬ 
porting  the  sewage  and  paper  mill  effluent  upstream  to  Harriman's  water 
intakes.  The  city  of  Harriman  was  eventually  forced  to  move  their  water 
intakes  further  upstream  in  the  Emory  River  arm  of  the  reservoir,  beyond 
the  influence  of  invading  density  currents. 

46.  Recent  research  on  reservoir  limnology  emphasizes  the  develop¬ 
ment  and  maintenance  of  longitudinal  gradients  in  reservoirs.  To  explain 
these  gradients,  Thornton  et  al.  (1981)  have  proposed  a  heuristic  model 
which  divides  a  reservoir  into  a  riverine  zone,  transition  zone,  and 
lacustrine  zone  (Figure  17).  The  definition  of  these  zones  is  depen¬ 
dent  on  how  reservoir  inflows  enter  and  move  through  reservoirs.  In 

the  riverine  zone,  current  velocities  are  decreasing  but  the  advective 
forces  are  still  sufficient  to  maintain  a  well-mixed  environment.  In 


Figure  17.  Riverine,  transitional,  and  lacrustrine  zones 
in  a  reservoir  (after  Thornton  et  al.  1981). 

the  transition  zone,  the  bouyancy  forces  begin  to  dominate  over  the 
advective  forces  and  the  inflow  plunges.  The  upstream  and  downstream 
boundaries  of  the  transition  zone  may  correspond  with  the  location  of 
the  plunge  point  under  low-flow  and  high-flow  conditions,  respectively; 
this  zone  is  also  the  zone  of  sedimentation.  In  the  lacustrine  zone, 
buoyancy  forces  dominate  and  inflows  move  through  the  reservoir  in  well- 
defined  horizontal  layers  as  interflows  and  underflows. 

Discussion 


47.  The  literature  review  indicates  that  an  understanding  of 
reservoir  density  currents  and  inflow  mixing  processes  is  essential  to 
understanding  reservoir  water  quality,  especially  with  respect  to  the 
Importance  of  storm  events.  The  review  also  indicates  that  the  analy¬ 
sis  of  density  currents  requires  an  in-depth  knowledge  of  inflow  charac¬ 
teristics,  ambient  reservoir  conditions,  and  reservoir  geometry.  It  is 
therefore  important  to  review  these  factors  prior  to  evaluating  predic¬ 
tive  techniques  for  density  currents. 

48.  The  important  attributes  of  inflow  are  flow  rate  and 
density.  Both  are  highly  variable.  Although  inflows  to  many  CE  reser¬ 
voirs  are  gaged  or  easily  back-calculated  from  outflows  and  changes  in 
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storage,  it  is  the  inflow  momentum,  rather  than  the  total  discharge,  that 
is  important  when  analyzing  inflow  density  currents.  The  inflow  momentum 
is  usually  difficult  to  assess  and  may  not  be  directly  related  to  the 
total  inflow  rate,  since  it  varies  as  the  width  and  thickness  of  the  flow 
zone  change.  The  inflow  density  is  a  function  of  temperature,  TDS,  and 
SS.  Inflow  temperatures  vary  with  hydrometeorological  conditions  and  are 
easily  measured  or  computed  from  heat  budget  considerations  (Ford  and 
Ford  1982) .  TDS  concentrations  can  increase  or  decrease  with  flow  and 
are  usually  related  to  specific  conductance,  which  is  easily  measured. 

SS  concentrations  can  sometimes  be  correlated  with  flow,  but  their  den¬ 
sity  contribution  is  complicated  by  settling  and  deposition  as  the  inflow 
enters  a  reservoir.  Although  the  inflow  characteristics  of  flow  rate 
and  density  are  both  measurable  and  predictable,  the  quantification  of 
inflow  momentum  and  SS  concentrations  is  not  a  simple  task. 

49.  The  ambient  reservoir  conditions  that  are  important  for  den¬ 
sity  current  analysis  include  density  stratification,  current  structure, 
and  turbulence  levels.  Ambient  reservoir  conditions  are  also  highly 
dynamic  but,  with  the  exception  of  temperature,  they  are  neither  easily 
measured  nor  predicted.  The  internal  current  structure  and  turbulence 
level  depend  on  such  ungovernable  factors  as  wind  speed,  variability,  and 
direction,  and  on  such  project  operation  factors  as  discharge  rate  and 
withdrawal  depth. 

50.  Irregular  reservoir  geometry  also  complicates  the  analysis  of 
density  currents  since  widths,  depths,  and  cross-sectional  areas  do  not 
vary  uniformly.  The  characteristic  dendritic  shape  of  reservoirs  usual] 
indicates  that  more  than  one  tributary  enters  the  reservoir  at  different 
locations  and  that  the  many  coves  and  embayments  may  be  isolated  and  no' 
part  of  the  zone  of  conveyance.  In  addition,  islands  may  produce  the 
bifurcation  of  an  inflow  density  current.  Other  factors  that  influence 
the  movement  of  density  currents  and  must  therefore  be  considered  are 
bottom  roughness  caused  by  uncut  timber,  for  example,  and  man-made 
obstacles  such  as  low-water  bridges,  earthen  embankments,  and 
cofferdams. 
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51.  In  summary,  when  evaluating  predictive  techniques  for  reser 
voir  density  currents,  one  must  realize  that  assumptions  such  as 
uniform  flow,  steady  state,  constant  width,  constant  density,  and  two 
layered  systems  are  not  representative  of  actual  conditions. 
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PART  III:  PREDICTIVE  TECHNIQUES 

52.  Numerous  simplified  techniques  have  been  proposed  in  the  lit¬ 
erature  to  predict  the  characteristics  of  density  currents.  These  are 
reviewed  and  evaluated  in  this  section  using  field  data  collected  from  a 
number  of  sources  including  Johnson  et  al.  1983  and  Kennedy  et  al.  1983. 
The  predictive  techniques  for  overflows,  plunge  point  location,  initial 
mixing  and  entrainment,  underflows,  and  intrusions  or  interflows  are 
discussed.  Use  of  the  equations  recommended  for  predicting  plunge  point 
location,  initial  mixing,  interflow  thickness,  and  propagation  speed  is 
demonstrated  in  Appendix  B,  employing  examples  from  DeGray  Lake,  Arkan¬ 
sas.  Appendix  C  outlines  an  inflow  algorithm  for  predicting  inflow 
placement,  plunge  point  depth,  initial  mixing,  underflow  entrainment, 
and  intrusion  thickness,  which  is  to  be  used  in  a  one-dimensional  water 
quality  model  such  as  CE-QUAL-R1  (Environmental  Laboratory  1982) . 

Overflows 

53.  Whenever  inflowing  water  is  less  dense  than  the  surface  water 
of  the  reservoir,  the  inflowing  water  will  float  on  the  water  surface 
and  an  overflow  occurs.  Several  problems  analogous  to  a  density  over¬ 
flow  have  been  studied  intensively.  These  include  thermal  discharges 
(e.g. ,  Adams  et  al.  1975,  Safaie  1979,  Harleman  1975)  and  oil  slicks 
(e.g.,  Wilkinson  1972,  1973).  In  addition,  Kao  (1977)  showed  that  a 
dense  front  penetrating  into  an  infinitely  deep  layer  of  lighter  fluid 
is  analogous  to  an  overflow  situation. 

54.  Three  factors  must  be  considered  when  analyzing  overflows; 
these  are:  the  separation  or  lift-off  point,  the  extent  of  horizontal 
spreading,  and  the  propagation  speed. 

Separation  point 

55.  The  separation  or  lift-off  point  is  that  point  where  the 
buoyancy  forces  exceed  the  advective  forces  and  the  inflowing  water 
floats  on  the  water  surface  (Figure  1).  The  separation  point  is 


33 


r 


analogous  to  the  plunge  point;  relationships  for  the  plunge  point  should 
also  be  valid  for  the  separation  point  (see  next  section  on  plunge 
points) . 

56.  In  a  series  of  laboratory  experiments  using  a  point  discharge 

into  a  wide  tank  over  a  sloping  bottom,  Safaie  (1979)  found  the  follow- 

* 

ing  empirical  equation  for  the  separation  depth  h  : 


h  =  0.914h.(F.) 
P  1  1 


where 

11^  =  separation  depth 
h^  =  inflow  depth 
and  the  inflow  Froude  number: 


where 


inflow  velocity 

density  of  receiving  water 

density  of  inflowing  water 


(pa  -  Pi> 


density  difference  between  receiving  water 
and  inflow 


g  *  gravitational  acceleration 

Equation  (1)  predicted  separation  depths  which  agreed  well  with  labora- 
2 

tory  experiments  (R  *  0.93)  and  was  verified  against  data  from  other 


Symbols  used  in  this  report  are  defined  once  where  used  in  the  main 
text  and  appendices  and,  for  convience,  are  listed  and  redefined  in 
Appendix  D.  Unless  specifically  stated,  all  equations  presented  in 
this  report  are  dimensionally  consistent,  and  either  metric  (SI)  or 
U.  S.  Customary  units  can  be  used  provided  the  user  is  consistent  in 
the  application. 
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investigations  and  field  studies  on  thermal  discharges.  Equation  (1)  is 
valid  for  F  >  1.2  .  If  F^  <  1.2  ,  then  h  <  h^  indicating  the  point 
of  separation  moves  up  into  the  inflowing  river  like  a  salt  wedge  in  an 
estuary.  This  result  is  consistent  with  the  work  of  Stigebrandt  (1978). 
Under  these  circumstances  (i.e.  F^  <  1.2  ),  the  separation  point  can  be 
assumed  to  occur  at  the  upper  end  of  the  reservoir. 

Horizontal  spreading 

57.  The  extent  of  horizontal  spreading  is  also  dependent  on  the 
inflow  Froude  number,  .  Safaie  (1979)  found  that  for  <  3, 
the  flow  was  mainly  an  unsteady  buoyancy-driven  spread  of  warm  water 
over  the  water  surface  in  all  directions.  Therefore,  for  F  <  3, 
it  can  be  assumed  that  an  overflow  will  be  well-mixed  laterally  for  all 
but  the  most  abrupt  entrance  geometry.  For  F  >  3  ,  the  initial 
spreading  was  controlled  by  the  jet  momentum  and  discharge  character¬ 
istics.  In  jets,  the  outer  boundaries,  regardless  of  their  definition, 
expand  proportionally  with  distance,  i.e., 

W  *  W  +  ax  (3) 

x  o 


where 

W  »  width  at  distance  x 
x 

W  *  initial  width 
o 

x  *  distance  from  source 

a  *  0.16  from  laboratory  experiments  for  plane  jets  (Fischer  et 
al.  1979) 

For  shallow  bodies  of  water,  bottom  interface  may  cause  greater  spread¬ 
ing  than  predicted  from  Equation  (3).  This  situation  will  probably  not 
occur  in  most  reservoirs. 

Propagation  speed 

58.  Kao  (1976,  1977)  showed  that  a  dense  front  penetrating  into  a 
infinitely  deep  layer  of  lighter  fluid  is  analogous  to  an  overflow 
situation;  since  this  problem  has  been  studied  extensively  it  will  be 
examined  first.  Von  Karman  (1940)  was  the  first  to  derive  an  expression 


for  the  steady  propagation  velocity  of  a  dense  front  into  a  quiescent 
lighter  fluid.  Assuming  a  perfect  fluid  and  irrotational  flow,  von 
Karman  used  Bernoulli's  equation  and  the  hydrostatic  relationship  to 
obtain 


u.  =  y2  ghs^  =  i.4iyghs^i. 


where 

ug  =  propagation  speed  of  the  front 

hg  =  thickness  of  the  underflowing  current 

59.  Benjamin  (1968)  showed  von  Karman' s  reasoning  was  incorrect 
because  a  breaking  headwave  (Figure  18)  was  necessary  to  balance  forces 
for  steady  propagation.  Without  the  breaking  headwave,  the  excess  flow 
force  per  unit  width. 


1  4ci  .2 

1—  8  hs 
a 


due  to  an  excess  in  hydrostatic  pressure  cannot  be  balanced.  A  drag 
force  (i.e.  form  drag)  due  to  the  velocity  deficit  in  the  wake  region 


STAGNATION 
POINT  y. 


behind  the  head  wave  must  be  considered.  Benjamin  rederived  the  propa¬ 
gation  speed  as  a  function  of  relative  depths  and  densities,  showing 
that  von  Karman's  result.  Equation  (4),  is  approached  when  the  ambient 
fluid  becomes  infinitely  deep. 

60.  The  speed  at  which  a  dense  front  propagates  into  a  lighter 
homogeneous  fluid  along  a  solid  boundary  has  been  studied  in  numerous 
laboratory  experiments  (Keulegan  1957,  Middleton  1966,  and  Wood  1966, 
among  others).  In  all  cases,  the  coefficient  C  in  the  equation 

u  =  C 
s 

was  found  to  be  less  than  the  theoretical  value  of  1.41  (See  Equation 
4).  In  a  review  of  these  experiments,  Benjamin  (1968)  found  that  the 
typical  value  for  C  was  approximately  1.1.  According  to  Benjamin, 
this  reduced  value  results  from  a  modified  piezometric  pressure  at  the 
stagnation  point,  which  in  practice  is  not  located  at  the  boundary  as 
assumed  with  inviscid  theory,  but  is  located  a  small  distance  up  from 
the  boundary  (Figure  18). 

61.  Laboratory  experiments  have  also  shown  that  stratification  in 
the  ambient  fluid  decreases  the  size  of  the  head  (Wilkinson  and  Wood  1972) 
which  explains  why  density  heads  are  not  always  observed  in  atmospheric 
flows  and  lakes.  In  addition,  Middleton  (1966)  found  that  the  shape  of 
the  front  remains  the  same  for  bottom  slopes  up  to  0.06.  This  discovery 
reinforces  the  fact  that  the  processes  controlling  the  shape  and  move¬ 
ment  of  the  front  are  distinct  and  different  from  the  frictional  forces 
that  control  the  uniform  flow  behind  the  head. 

62.  Following  Benjamin  (1968),  Kao  (1977)  derived  a  general  ex¬ 
pression  for  the  steady  propagation  speed  of  a  nonuniform  density  cur¬ 
rent  into  a  stratified  ambient  environment.  Kao  showed  that  (a)  the 
effect  of  the  ambient  stratification  was  to  increase  the  speed  of  the 
density  current  and  (b)  the  speed  of  warmwater  surface  current  was  in¬ 
dependent  of  mixing  and  depended  only  on  the  inflow  Froude  number;  i.e.. 


us  2  1/3 

^  '  <2/Fi> 


(6) 


63.  Koh  (1976)  analyzed  the  time-dependent  spreading  of  a  buoyant 

fluid  on  the  surface  of  a  heavier  fluid.  He  assumed  the  shape  of  the 
spreading  fluid  was  geometrically  similar  from  one  instant  to  the  next, 
and  divided  the  problem  into  three  time  scales:  initially,  there  is  a 
transient  start-up  period;  then  there  is  an  intermediate  time  where  the 
inertial  forces  balance  the  buoyancy  driving  forces;  and  finally,  for 
large  times,  interfacial  shear  is  dominant.  For  practical  problems, 

the  intermediate  time  scale  is  of  interest,  and  Koh  determined  the  fol¬ 
lowing  time  (t)  dependence  for  overflow  length  and  thickness: 

length 

2/3 

two-dimensional  instantaneous  release  t 

two-dimensional  continuous  release  t 

1/2 

axisymmetric  instantaneous  release  t 

3/4 

axisymmetric  continuous  release  t 

Therefore,  for  the  two-dimensional  continuous  release  situation  the 
overflow  length  is  proportional  to  time  t  and  the  thickness  is  con¬ 
stant.  For  large  times,  however,  Koh  showed  that  the  thickness  is 
proportional  to  t^"*  and  would  therefore  increase  indefinitely.  In 
contrast,  the  continuous  release  axisymmetric  case  tends  to  a  constant 
thickness  for  large  times.  The  third  dimension  is  sufficient  to  pre¬ 
vent  blocking  of  the  flow. 

64.  Assuming  shear  is  small  and  the  flow  regime  is  governed  by  a 
balance  of  inertial  and  buoyant  forces,  Koh  (1976)  solved  the  two- 
dimensional  continuous  release  problem  and  obtained: 


thickness 


u 

s 


(7) 
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where 


L  =  overflow  length 
s 

q  =  unit  width  discharge 

C  =  coefficient  to  account  for  nonhydrostatic  pressure  distribu- 
**  tion 

Using  C  =0.5  from  laboratory  experiments  for  a  surface  source, 
Equation  (7)  reduces  to 


u  =  1.03 
s 


(8) 


and  assuming  uniform  flow  (q^  =  ushg) 


u 

s 


(9) 


Since 


q .  t  =  -r  L  h 
l  4  s  s 


for  an  overflow  suaped  like  one  quarter  of  an  ellipse  (Koh  1976) ,  then 


h  =1.24 
s 


qi 


Ap. 

l 


1/3 


(10) 


The  coefficient  in  Equation  (9)  (1.04)  is  similar  to  the  value  of  1.1 
recommended  by  Benjamin  (1968)  (Equation  5). 

Summary 

65.  Several  factors  complicate  the  analysis  of  overflows.  First, 
three-dimensional  effects  must  be  considered.  Since  the  motive  force 
in  an  overflow  is  the  excess  hydrostatic  pressure,  the  current  will 
spread  in  any  direction  not  obstructed  by  boundaries.  Second,  an  over¬ 
flow  will  usually  result  from  density  differences  due  to  temperature. 


39 


This  temperature  difference  is  quickly  eliminated  by  heat  transfer  at 
the  air-water  interface.  Third,  any  wind  shear  can  (a)  direct  the  over¬ 
flow  into  a  cove  or  prevent  it  from  moving  downstream  and  (b)  enhance 
horizontal  dispersion.  Fourth,  vertical  mixing  resulting  from  wind 
shear  and/or  convective  cooling  can  distribute  the  current  throughout 
the  water  column. 

66.  The  overflow  speed  and  thickness  can  be  estimated  from  Equa¬ 
tions  (9)  and  (10),  respectively. 

Plunge  Point  Analysis 

67.  If  the  density  of  the  inflowing  river  water  is  greater  than 
the  surface  water  of  the  reservoir,  the  inflow  will  plunge  beneath  the 
water  surface.  The  location  of  the  plunge  point  is  determined  by  a 
balance  of  the  stream  momentum  (i.e.,  dynamic  force);  the  pressure 
gradient  across  the  interface  separating  the  river  water  from  the 
reservoir  water  (i.e.,  static  force);  and  the  resisting  shear  forces  at 
the  bed,  interface,  and  water  surface.  Since  flow  changes  in  the 
vicinity  of  the  plunge  point  occur  over  distances  small  compared  to 
normal  cross-section  changes,  the  flow  is  rapidly  varying  and  therefore 
independent  of  slope. 

68.  Singh  and  Shah  (1971)  provided  an  excellent  description  of 
the  flow  in  the  vicinity  of  a  plunge  point  (Figure  19).  The  location 
maximum  velocity  gradually  changes  from  near  the  water  surface  in  a 


Figure  19.  Flow  in  the  vicinity  of  the  plunge  zone 
(after  Singh  and  Shah  1971). 
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riverine  zone  to  near  the  bottom  in  the  reservoir  backwater  as  the 
static  force  begins  to  dominate.  At  the  plunge  point,  the  flow  is 
concentrated  near  the  bed. 

69.  For  small  slopes  and  Froude  numbers,  under  more  or  less  steady 
state,  the  sharp  line  usually  delineating  the  plunge  point  indicates 
that  little  mixing  is  occurring.  However,  for  larger  Froude  numbers  and 
unsteady  conditions,  mixing  is  probably  significant  and  proportional  to 
the  energy  loss  that  occurs  in  the  vicinity  of  the  plunge  point. 

70.  Various  investigators  have  proposed  simple  models  for  the 
plunge  location  (Singh  and  Shah  1971,  Wunderlich  and  Elder  1973,  Savage 
and  Brimberg  1975,  Hebbert  et  al.  1979,  Jain  1981,  Akiyama  and  Stefan 
1981).  Some  of  these  have  been  based  on  experimental  work  (Singh  and 
Shah  1971)  and  others  on  theoretical  work  (Hebbert  et  al.  1979,  Savage 
and  Brimberg  1975,  Jain  1981,  Akiyama  and  Stefan  1981).  Assumptions 
vary  with  the  investigator,  but  typically  include  steady  state,  constant 
width,  no  mixing,  uniform  velocity  profiles,  and  constant  reservoir 
density . 

Singh  and  Shah 

71.  Singh  and  Shah  (1971)  conducted  an  experimental  study  of  the 
plunging  phenomenon  using  a  tilting  flume  with  salt  water  flowing  into  a 
reservoir  filled  with  tap  water.  Using  the  method  of  least  squares, 
Singh  and  Shah  related  the  plunge  point  depth,  h^  ,  to  the  critical 
depth  by  the  equation: 

h  -  0.0185  +  1.3  l  — —  I  (11) 

P  \  8£i  / 

where 

hp  ■  plunge  point  depth,  m 

■  6p^/pa  “  relative  density  difference 
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In  this  study,  the  unit  discharge  was  varied  from  0.5  x  10  to  135 
-4  2 

x  10  m  /s,  the  flume  slope  was  varied  from  0.005  to  0.02,  and  the  den¬ 
sity  difference  was  varied  from  0.5  to  13.0  g/i,.  The  Reynolds  number 
ranged  from  600  to  11,000.  For  practical  applications  where  plunge 
point  depths  are  the  order  of  meters,  the  constant  (0.0185),  can  be 
assumed  small  and  Equation  (11)  becomes: 


h  = 
P 


(12) 


72.  Singh  and  Shah  also  investigated  the  plunging  phenomenon  analy¬ 
tically.  Applying  the  momentum  principal  across  the  transition  region 
(Figure  20),  they  were  able  to  derive  an  equation  similar  to  the  equa¬ 
tion  for  conjugate  depths  of  a  submerged  hydraulic  jump.  Substituting 
experimental  values  for  friction,  they  reduced  the  momentum  equation  to: 


Figure  20.  Schematic  of  plunge  zone. 


3.1 


(hj_  +  h2)  (h1h2) 


where 

■  observed  plunge  depth 
h2  *  thickness  of  underflow 

Assuming  h^  is  similar  to  h2  (h^  'v  h2  ^  h^) ,  then  this  equation 
reduces  to: 


(A 

\  -  1.16^  J  (13) 


which  is  similar  to  Equation  (12). 

Wunderlich  and  Elder 

73.  Following  Harleman  (1961),  Wunderlich  and  Elder  (1973)  assumed 
that  conditions  at  the  plunge  point  were  governed  by  a  densimetric  Froude 
number 


F 

P 


A  .  „  .  ,1/2 

(ge^p) 


(14) 


where 


*  discharge 

A  *  h  W  ■  cross-sectional  area 
P 

W  *  width 

74.  Using  data  from  the  Natahala  River  inflow  to  Fontana  Reser¬ 
voir,  Wunderlich  and  Elder  found  F^  =  0.5  at  the  plunge  point.  Sub¬ 
stituting  Fp  *  0.5  into  Equation  (14)  and  rearranging  terms  yields 
the  plunge  point  depth: 


h 


P 


4.0  <VA)' 
8e< 


(15) 
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Savage  and  Brimberg 

75.  Savage  and  Brimberg  (1975)  analyzed  plunging  phenomena  in  two 
ways:  the  first  was  based  on  conservation  of  energy  and  the  second  on 
gradually  varied  flow  theory  in  a  two-layered  stratified  system. 

76.  In  the  simple  energy  balance,  interfacial  and  bed  friction  as 
well  as  slope  were  neglected.  The  flow  in  the  upper  layer  was  assumed 
to  be  zero.  The  Froude  number  at  the  plunge  point  was  found  to  be  0.5, 
which  resulted  in  a  plunge  point  depth  of: 

2  v  1/3 


h  =1.59 
P 


(-t) 


(16) 


This  equation  was  compared  with  the  experimental  data  of  Singh  and  Shah 
(1977)  and  was  found  to  fit  the  data  set  as  well  as  the  empirical  rela¬ 
tionship  (Equation  11). 

77.  In  the  two-layered  analysis,  Savage  and  Brimberg  used  the  one¬ 
dimensional  equations  of  motion  for  gradually  varying  flow  presented  by 
Schijf  and  Schonfeld  (1953).  These  equations  included  bed  slope  and 
interfacial  and  bed  friction  and  were  used  to  define  the  shape  of  the 
interface.  For  the  case  where  flow  in  the  lower  layer  reaches  normal 
depth  downstream  of  the  plunge  point  (i.e.,  an  M2-type  profile),  the 
equation  was  integrated  numerically  to  determine  the  Froude  number  as  a 
function  of  slope  and  roughness.  That  is. 


F 

P 


2.05  (  S  \  0,478 
d+a)  \  ffa  / 


(17) 


where 


a  *  f./f,  =  ratio  of  interfacial  to  bed  friction 
1  b 

f.  =  bed  friction  coefficient 
b 

f^  *  interfacial  friction  coefficient 
S  ■  bed  slope 

Savage  and  Brimberg  give  ranges  of  0.2  to  0.8  for  a  and  0.01  to  0.09 
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for  f, 
o 

tion 


The  value  for  F  can  then  be  substituted  into  the  equa- 
P 


to  determine  the  plunge  point  depth  as  a  function  of  slope  and  friction 
as  well  as  discharge  and  density  difference. 

78.  Since  Singh  and  Shah  did  not  present  densimetric  Froude  num¬ 
bers  for  their  experiments.  Savage  and  Brimberg  estimated  the  Froude 
numbers  to  be  between  0.3  and  0.8.  Substitution  of  these  limits  into 
Equation  (17)  banded  the  data  of  Singh  and  Shah  and  accounted  for  their 
previous  scatter. 

Hebbert  et  al. 

79.  Instead  of  considering  a  rectangular  cross  section  of  unit 
width,  Hebbert  et  al.  (1979)  considered  a  triangular  cross  section  with 
half  angle  $  (Figure  21).  Starting  with  the  equations  for  conserva¬ 
tion  of  volume  and  momentum,  Hebbert  et  al.  related  the  downstream 

normal  densimetric  Froude  number  F  to  the  ratio  between  the  down- 

n 

stream  normal  depth  of  flow  h  and  the  plunge  point  depth  h 

n  p 

(i.e.,  C  *  hR/hp  ).  The  plunge  point  depth  can  then  be  calculated, 

knowing  the  densimetric  Froude  number  and  downstream  normal  depth  of 

flow.  Hebbert  et  al.  related  F  to  reservoir  characteristics: 

n 


n 


sin  S) 


(19) 


where  CQ  «*  (f^  +  f^)/4  =  drag  coefficient. 

Equation  (19)  states  that  the  downstream  normal  densimetric  Froude 
number  is  a  function  of  reservoir  characteristics  only  and  therefore 
a  constant  for  a  given  reservoir.  The  downstream  normal  depth  of  flow 
is  then  computed  from 


80.  If,  instead  of  a  triangular  cross  section,  Hebbert  et  al.  had 

assumed  a  rectangular  cross  section  of  unit  width,  then  the  difference 

between  the  plunge  point  depth  and  downstream  normal  depth  of  flow  is 

smaller  (i.e.,  C  =  ^n  =  0.99). 

h 

P 

Jain 

81.  Jain  (1981)  critically  examined  the  analysis  of  Savage  and 
Brimberg  (1975) .  He  showed  that  the  simple  energy  conserving  model  was 
physically  unrealistic  since  energy  dissipation  near  the  plunge  point 
was  required  for  the  flow  mechanism,  a  demonstration  which  also  invali¬ 
dates  the  assumption  that  the  plunge  point  is  a  stagnation  point. 

82.  Jain  also  examined  the  gradually  varied  two-layer  flow  analy¬ 
sis  of  Savage  and  Brimberg  (1975)  by  starting  with  the  equations  of 
Schijf  and  Schonfeld  (1953)  and  assuming  the  mean  velocity  in  the  upper 
layer  is  zero.  He  concluded  that  the  interfacial  profile  calculation 
was  sensitive  to  the  direction  of  numerical  integration.  The  equations 
should  be  integrated  in  the  direction  in  which  the  hydraulic  control  is 
acting.  For  mild  slopes,  which  are  characteristic  of  most  reservoir 
systems,  the  hydraulic  control  is  downstream  and  the  integration  should, 
therefore,  start  downstream,  not  at  the  plunge  point. 

83.  Based  on  the  gradually  varied  two-layer  flow  analysis,  Jain 
proposed  a  nondimensional  formula  for  the  plunge  point  depth 


N 

P 


,  0.126  „  0.024 

1.  DA  N 

C 


(23) 


where 


N 

P 


h 

n 


h  /h 
P  n 


normal  depth 


X 

N 


«/(l  + 
h  /h 


c  n 


a) 


47 


=  critical  depth 


(A) 


Equation  (23)  can  be  rearranged  to  solve  for  the  plunge  depth  directly 


0.325 


h  =  0.814 
P 


Akiyama  and  Stefan 

84.  Akiyama  and  Stefan  (1981)  proposed  a  conceptual  flow  model  for 
the  plunge  point  which  depends  on  inflow  conditions  (densimetric  Froude 
number),  downstream  conditions  (channel  slope,  roughness,  and  width),  and 
mixing.  Starting  with  the  integrated  momentum  equation  and  assuming 
steady  flow,  they  determined  the  plunge  depth  for  a  mild  slope: 


h  =  \ 
P  2 


where 


K  =  ~  constant 

f  =  f,+  f.  =  total  friction  factor 
t  b  l 


S2  e. 


k 


Pa  -  P(z) 

g  -  dz  =  density  deficiency 


ek  =  (°k_  Pj/ Pa  =  Ap  /p 
b  baa  b  a 

=  underflow  density 

h^  =  depth  of  underflow  downstream  of  plunge  point 
Y  »  initial  mixing  coefficient 
p(z)  *  density  variation  with  depth,  z 
and  for  a  steep  slope 
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I  (1*  si*y' 


(1  +  sA)' 


(1  +  Y>‘ 


■m'%) 


1/3 


(1  +  Y) 


(26) 


where 


z  dz  =  density  deficiency 


Since  most  reservoirs  are  characterized  by  a  mild  slope  (i.e.,  S 

_3 

10  ),  only  the  solution  for  mild  slope  will  be  discussed  hereafter. 

85.  The  equation  for  plunging  depth  (Equation  25)  on  a  mild  slope 
consists  of  three  parts.  The  first  term  on  the  right  side  of  Equation 
(25)  corresponds  to  depth  ratio  (tWh^)  .  For  mild  slopes,  the  ratio 
(S/K)  is  small  compared  to  1,  and  Equation  (25)  reduces  to 


•«*’>  (f)'/3(Al) 


1/3 


(27) 


Equation  (27)  indicates  that  when  mixing  is  zero  (y  =  0)  ,  the  plunge 
depth  is  the  normal  depth  downstream  of  the  plunge  point.  Mixing 
(i.e.,  y  )  causes  the  plunge  depth  to  be  larger  than  the  normal  depth. 
Equation  (27)  can  be  further  reduced  to 


(28) 


if  it  is  assumed  that  =  0.75  (Ellison  and  Turner  1959)  since 
K  -  ft/S2  . 

Discussion 

86.  All  of  the  models  locate  the  plunge  point  by  calculating  the 
hydraulic  depth  at  which  the  inflow  plunges,  which  can  be  put  in  the 
general  form: 
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r 


h 

P 


Wcs£i 


1/3 


(29) 


where 

h  =  hydraulic  depth  at  the  plunge  point 

Q  =  original  inflow  rate 

W  =  width  of  the  zone  of  conveyance  (the  zone  where  the  density 
C  current  flows) 

g  =  acceleration  due  to  gravity 

e.  =  relative  density  difference  between  the  lake  and  inflow  waters 

F  =  critical  densimetric  Froude  number 
P 

Akiyama  and  Stefan  (1981)  included  a  mixing  term.  The  formulation 
of  Hebbert  et  al.  (1979)  differs  slightly  from  Equation  (29)  because 
they  assumed  a  triangular  instead  of  a  rectangular  cross  section. 

87.  Before  discussing  the  use  of  Equation  (29),  two  important 

points  need  to  be  emphasized.  First,  Equation  (29)  solves  for  the 

hydraulic  depth  at  the  plunge  point  (i.e.,  h^  =  cross-sectional  area  * 

reservoir  width),  not  the  average  or  maximum  depth.  The  hydraulic  depth 

is  not  readily  available  from  bathymetric  maps;  it  must  be  calculated 

from  reservoir  cross-section  data.  The  triangular  cross  section  used 

in  the  Hebbert  et  al.  (1979)  formulation  provides  a  unique  relationship 

between  the  reservoir  width,  hydraulic  depth,  and  distance  along  the 

longitudinal  axis.  Use  of  this  method  therefore  has  distinct  advantages 

over  the  rectangular  shaped  cross  section.  Second,  Equation  (29)  was 

developed  for  constant-width  channels  with  the  flow  rate  expressed  as 

discharge  per  unit  width  (i.e.,  q^  =  Q^/W  ).  Since  the  width  is  not 

constant  for  most  reservoirs,  Equation  (29)  cannot  be  solved  directly 

for  h  .  In  addition,  the  width  W  is  not  the  total  width  but  the 
P 

width  of  conveyanc  .  zone.  As  stated  earlier,  DeGray  Lake  is  not  charac¬ 
terized  by  complete  lateral  mixing.  Even  in  more  uniformly  shaped 
reservoirs,  density  flows  tend  to  follow  the  old  river  channel 
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(Kennedy  et  al.  1981). 

88.  Values  of  F  typically  vary  from  approximately  0.1  to  0.7 

for  a  mild  slope.  Savage  and  Brimberg  (1975)  found  F  =  0.5  based  on  a 

P 

simple  energy  balance.  Wunderlich  and  Elder  (1973)  also  found  that  this 

value  correctly  predicted  plunge  depths  in  Fontana  Reservoir,  yet  there 

is  is  no  reason  to  believe  that  F  =  0.5  should  be  a  universal  con- 

P 

stant.  Using  a  momentum  balance,  Savage  and  Brimberg  (1975)  and  Akiyama 
and  Stefan  (1981)  formulated  F^  in  terms  of  bed  slope,  bed  roughness, 
and  interfacial  roughness.  Once  F^  is  determined,  the  plunge  depth 
can  be  calculated  from  Equation  (29)  in  an  iterative  fashion  by  assuming 
a  width  and  calculating  a  depth.  The  widths  are  continually  updated 
until  they  correspond  with  an  actual  width  occurring  in  the  reservoir 
at  the  calculated  hydraulic  depth  for  the  plunge  point.  The  major 
problem  in  using  Equation  (29)  is  that  the  width  is  not  always  uniquely 
related  to  the  hydraulic  depth,  especially  if  plug  flow  does  not  occur. 

89.  Using  the  data  from  the  storm  events  and  dye  studies  sum¬ 
marized  in  Table  1,  plunge  point  depths  were  calculated  with  the  equa¬ 
tions  recommended  by  Singh  and  Shah  (1971),  Wunderlich  and  Elder  (1973), 
Savage  and  Brimberg  (1975),  Hebbert  et  al.  (1979),  Jain  (1981),  and 
Akiyama  and  Stefan  (1981)  (i.e..  Equations  12,  15,  18,  22,  24,  and  28, 
respectively).  Site-specific  parameters  were: 


DeGray  Lake 


a 

Y 


West  Point  Reservoir 
-4 


1.2  x  10"3 

8.0  x  10 

0.004 

0.004 

0.02 

0.02 

0.024 

0.024 

0.2 

0.2 

0 

0 

O 

m 

oo 

87° 

Excluding  the  results  from  Singh  and  Shah  and  Akiyama  and  Stefan,  the 
plunge  point  predictions  (Table  2)  were  generally  within  +  1  m  or  20 


51 


Summary  of  Inflow  Studies 
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Specific  conductance 


percent,  provided  the  width  of  the  zone  of  conveyance  was  used.  The 
predictions  listed  in  Table  2  for  Hebbert  et  al.  are  one-half  the  result 
derived  from  Equation  (22),  since  Equation  (22)  predicts  the  maximum 
plunge  depth  not  the  hydraulic  depth  (i.e.,  for  a  triangular  cross 
section,  the  hydraulic  depth  is  one-half  the  maximum  depth).  The  pre¬ 
dictions  from  the  Singh  and  Shah  equation  tended  to  be  too  small  while 
the  predictions  from  the  Akiyama  and  Stefan  equation  tended  to  be  too 
large.  If  the  total  friction  coefficient  f  was  reduced  to  0.006, 
then  the  Akiyama  and  Stefan  predictions  were  almost  identical  to  the 
predictions  from  the  Savage  and  Brimberg  equation.  This  finding  is  not 
unexpected  since  the  reduced  forms  of  Equations  (18)  and  (28)  are 
similar. 

90.  In  general,  the  plunge  point  depth  predictions  were  better  for 
the  steady  base  flow  conditions  than  for  the  unsteady  storm  events. 
Plunge  depths  for  storm  events  can  be  computed,  however,  provided  the 
maximum  flow  and  inflow  density  at  the  time  of  maximum  flow  were  used. 
Attempts  to  use  other  conditions  such  as  mean  daily  values  or  averaged  , 
storm  values  were  unsuccessful. 

91.  Sensitivity  analysis  on  the  six  predictive  equations  indicated 
the  results  were  more  sensitive  to  perturbations  in  flow  rate  Q 

and  conveyance  width  than  to  pertubations  in  density  differences 

e^  .  In  some  storm  events,  a  change  of  one  order  in  magnitude  of 

resulted  in  a  change  of  hydraulic  depth  prediction  of  less  than  1  m, 
while  a  50  percent  change  in  flow  (i.e.,  one-half  the  maximum  flow) 
resulted  in  a  change  of  2  m  or  more.  This  sensitivity  results  from 

Q.  and  W  being  squared  in  Equation  (29). 

I  c 

92.  It  is  recommended  that  the  equations  of  Savage  and  Brimberg 
(Equation  18),  Hebbert  et  al.  (Equation  22),  Jain  (Equation  24),  and 
Akiyama  and  Stefan  (Equation  28)  be  used  for  plunge  point  predictions 
since  they  consider  slope  and  friction  as  well  as  flow  and  density 
differences.  They  can  therefore  be  calibrated  to  specific  systems.  It 
is  also  recommended  that  all  four  equations  be  used  and  the  predictions 
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compared  when  a  new  or  different  reservoir  is  being  evaluated.  It  is 
possible  that  for  any  specific  reservoir  a  different  equation  may  be 
more  appropriate.  For  reservoirs  that  gradually  expand,  the  Hebbert  et 
al.  equation  is  recommended  since  there  is  a  unique  relationship  between 
width  and  longitudinal  distance  and  the  predicted  plunge  depths  corres¬ 
pond  to  the  maximum  reservoir  depths.  In  the  other  three  equations,  the 
hydraulic  plunge  depth  is  predicted  and  the  relationship  between 
hydraulic  depth  and  width  must  be  known.  It  must  be  noted,  however, 
that  the  hydraulic  depth  may  not  be  uniquely  related  to  longitudinal 
distance. 


Initial  Mixing  and  Entrainment 

93.  Initial  mixing  and  entrainment  are  two  types  of  mixing  phe¬ 
nomena  important  when  analyzing  density  currents.  Initial  mixing 
includes  the  cumulative  effects  of  all  mixing  processes  acting  in  the 
vicinity  of  the  plunge  point.  Entrainment  occurs  at  the  interface 
between  the  underflow  and  reservoir  water  after  the  inflow  has  plunged. 
Turbulence  generated  by  bottom  roughness  entrains  reservoir  water  into 
the  underflow. 

Initial  Mixing 

94.  At  the  plunge  point,  the  flow  changes  from  a  homogeneous  open- 
channel  flow  to  a  stratified  underflow.  Water  tends  to  pool  at  the 
plunge  point,  since  it  flows  into  the  plunge  zone  faster  than  it  flows 
out  creating  a  large  eddy  (Knapp  1942,  Elder  and  Wunderlich  1972,  Ford 
et  al.  1980).  Changes  in  flow  depth  and  direction  all  contribute  to 
mixing.  Since  the  mechanism  is  complex,  the  initial  mixing  is  usually 
incorporated  into  an  empirical  coefficient  y  defined  by 

Qbo  =  (1  +  (30) 


Q 


bo 


where 


revised  flow  rate  including  entrainment  downstream  of  the 
plunge  point 


The  density  of  the  inflow  is  decreased  by 


(31) 


where 

p,  =  revised  inflow  density  after  entrainment  at  the  plunge  point 
bo 

95.  Field  data  have  been  collected  by  Chasse  and  Slotta  (1972), 

Elder  and  Wunderlich  (1972),  and  Ford  et  al.  (1980).  Since  these  data 
may  also  include  the  effects  of  entrainment  into  the  underflow  (to  be 
discussed  in  the  next  section) ,  the  mixing  coefficient  will  be  called 
the  entrance  mixing  ratio  r  .  The  mixing  ratio  r  is  defined  by 
Equations  (30)  and  (31)  with  y=r ,  but  Q  and  p^  are  defined  far 
downstream  of  the  plunge  point.  In  some  instances,  they  are  the  flow 

rate  and  density  at  the  point  where  the  underflow  leaves  the  old  river 

channel  to  become  an  interflow. 

96.  Since  y  (or  r  )  is  expected  to  be  a  function  of  the  den- 
simetric  Froude  number,  mixing  ratios  from  Chasse  and  Slotta  (1972), 
Elder  and  Wunderlich  (1972),  and  Ford  et  al.  (1980)  are  plotted  versus 
densimetric  Froude  number  in  Figure  22.  Three  observations  are  worth 
noting: 

a.  The  data  are  widely  scattered, with  F  varying  from  0.1  to 

0.7  and  r  varying  from  0  to  0.8.  This  scatter  is  not  un¬ 

expected  considering  the  difficulties  encountered  when  col¬ 
lecting  field  data  of  this  nature  and  the  data  manipulations 
required  to  obtain  F  and  r  . 

b.  Mixing  ratios  for  specific  reservoirs  are  more  consistent 

and  appear  to  be  independent  of  F  .  For  example,  the  average 
mixing  ratios  for  DeGray  Lake,  Arkansas,  Chilhowee  Reservoir, 
Tennessee,  Norris  Reservoir,  Tennessee,  and  Hills  Creek  Reser¬ 
voir,  Oregon, are  0.12,  0.03,  0.38,  and  0.54,  respectively. 

£.  The  average  mixing  ratio  for  all  data  is  0.26  which  is  con¬ 
sistent  with  the  initial  mixing  coefficient  of  0.3  recommended 
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and  inflow  densimetric  Froude 


by  Akiyama  and  Stefan  (1981) . 

97.  Ryan  and  Harleman  (1971)  reported  mixing  ratios  obtained 
from  laboratory  experiments  to  vary  from  0.5  to  2.0.  Modification  of 
the  entrance  to  the  laboratory  flume  reduced  the  mixing  ratios  to 
between  0.1  to  0.5  indicating  that  the  geometry  of  the  inflowing  river 
may  be  an  important  factor. 

98.  It  is  apparent  that  little  is  known  about  inflow  mixing.  The 
field  data  indicate  that  a  constant  entrance  mixing  ratio  may  be  ap¬ 
propriate  for  specific  reservoirs,  but  physical  reasoning  would  indicate 
that  the  mixing  ratio  should  be  dependent  on  flow  and  the  densimetric 
Froude  number.  The  mixing  ratio  should  increase  as  the  advective 
forces  dominate  over  the  buoyant  forces.  A  relationship  between  y 

and  F  such  that  ^  increases  as  F  increases  is  the  empirical 
equation 

Y  =  1.2  F  -  0.2  (32) 

recommended  by  Jirka  and  Watanabe  (1980)  for  cooling  ponds  which 
is  valid  for  F  ±  0.167. 

Entrainment 

99.  After  the  inflow  plunges,  it  generally  follows  the  old  river 
channel  (thalweg)  as  an  underflow.  Turbulence  generated  by  bottom 
roughness  entrains  reservoir  water  into  the  underflow.  This  second  type 
of  mixing  has  been  analyzed  by  Ellison  and  Turner  (1959),  Hebbert  et  al . 
(1979),  and  Akiyama  and  Stefan  (1981).  Entrainment  rates  for  an  under¬ 
flow  can  be  calculated  from  field  data  using  conservation  of  volume: 


M  +  f-x(ubA)  =  EubB 


where 

ub  ■*  mean  velocity  of  the  underflow 
B  =  top  width  of  the  underflow 
A  ■  cross-sectional  area  of  flow 
E  “  the  entrainment  coefficient 


(33) 
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Equation  (33)  assumes  that  the  entrainment  coefficient  is  proportional 

to  the  underflow  velocity.  Elder  and  Wunderlich  (1972)  and  Hebbert  et 

-4 

al.  (1979)  calculated  entrainment  rates  on  the  order  of  10  .  Entrain- 

_3 

ment  rates  for  DeGray  Lake  were  determined  to  be  on  the  order  of  10 
The  larger  rate  found  in  DeGray  Lake  may  be  attributed  to  wind-induced 
mixing  at  the  bottom  of  the  epilimnion,  or  to  increased  roughness  caused 
by  trees  that  were  not  cleared  from  the  reservoir  prior  to  filling  (Ford 
and  Johnson  1981) . 

100.  The  entrainment  coefficient  is  normally  assumed  to  be  a 
function  of  the  Richardson  number  (or  densimetric  Froude  number). 

Using  experimental  data  from  Ellison  and  Turner  (1959),  Lofquist  (1960), 
Kato  and  Phillips  (1969),  and  their  own  experiments,  Ashida  and  Egashira 
(1977)  found 

E  *  0.0015  Ri-1  (34) 


where 

Ri 


where 


Apb 

hb 


*  density  difference  of  underflow  and  ambient 
=  depth  (thickness)  of  underflow 


Imberger  and  Patterson  (1981)  suggested 


E  = 


3  3/2p2 

n  CKCb  % 


(35) 


where 

3 

n  =  mixing  efficiency 

=  bottom  drag  coefficient 

3 

with  the  efficiency  n  =  3.2  based  on  the  field  studies  of  Hebbert 
et  al.  (1979)  and  Elder  and  Wunderlich  (1972).  This  mixing  efficiency 
is  a  factor  of  10  larger  than  the  mixing  efficiency  at  the  bottom  of  the 
mixed  layer  (Imberger  and  Patterson  1981).  Equation  (35)  has  the 
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advantage  of  including  the  bottom  drag  coefficient.  Assuming  = 
0.015  (Imberger  and  Patterson  1981),  Equation  (35)  becomes 

E  =  2.9  x  10”3  F3  =  2.9  x  10-3Ri_1 
b 

which  is  similar  to  Equation  (34). 

Underflows 


101.  After  the  inflowing  river  passes  through  the  plunge  point 
zone,  it  moves  down  the  submerged  river  channel  as  an  underflow.  The 
depth  (thickness)  of  the  underflow  is  determined  by  the  discharge,  the 
density  of  flow,  and  the  bed  and  interfacial  friction.  Gradually  vary¬ 
ing  flow  theory  can  be  used  to  describe  the  flow.  If  the  bottom  slope 
is  mild  (S  <  for  f t  =  0.02),  the  underflow  will  remain  subcritical 
and  uniform  flow  will  be  achieved.  If  the  bottom  slope  is  steep,  the 
flow  will  become  supercritical  downstream  of  the  plunge  point.  Since 
most  reservoir  bottom  slopes  are  mild,  only  mild  slopes  will  be  dis¬ 
cussed  herein.  The  steep  slope  or  supercritical  condition  is  discussed 
in  Akiyama  and  Stefan  (1981). 

102.  Following  Ellison  and  Turner  (1959)  ,  the  integral  method  for 
analyzing  the  two-dimensional  underflow  is  used  here  to  avoid  making 
assumptions  concerning  the  shape  of  the  velocity  and  density  profiles. 
If  hydrostatic  equilibrium  and  uniform  flow  are  assumed  and  if  the 
coordinate  system  is  taken  parallel  to  the  bottom  slope  such  that  the 
velocity  normal  to  the  boundary  is  assumed  small,  then,  for  an  elemental 
volume  of  unit  width,  conservation  of  momemtum  is 

d  (u3  h,  )  =  -f  -Id  (S  g'h3  cos  <0  +  S.g’h.sin  <*  (36) 

~r~  b  b  t  b  yr  —  lb  2  b 

dx  2  dx 

conservation  of  volume  is 


d_  (ubhb)  *  Eub  (37) 

dx 
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and  conservation  of  mass  is 


d_  (APbhbub)  =  0  (38) 

dx  * 


where 

=  depth  (thickness)  of  underflow 
f  =  total  drag  coefficient 
E  =  entrainment  coefficient 
0  =  bed  slope 

°b  -  pa 

g'  =  relative  acceleration  of  gravity  =  -  g 


S^,  “  profile  constants 

In  Equation  (36),  the  first  term  on  the  right  side  is  the  frictional 
drag,  the  second  term  is  the  difference  of  the  integrated  hydrostatic 
pressure,  and  the  last  term  is  the  gravitational  force  accelerating  the 
layer.  Following  Akiyama  and  Stefan  (1981),  Equations  (36),  (37)  and 
(38)  can  be  manipulated  to  obtain 


■  «2  -  kR1)E  -  S2R1  tan 


+  ft)(l  -  S1Ri)"1Ri'1  (39) 


b  dRi 


3  Ri  dx 


((1  +  ^RiJE  -  S2 


Ri 


tan  0  +  f  ) (1 


S1Ri) ~1Ri~^ 


(AO) 


where 


Ri 


g'h^  cos  0 


2 

“b 


(41) 


and  it  is  assumed  that  ,  S,.  ,  and  @  are  independent  of  distance. 

Equations  (39)  and  (40)  can  be  further  simplified  to 
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(42) 


r 


^b_  Jlb_dRi_]^  §  \  dRi 

dx  "  3Ri  dx  Ri  n.n-l  3  dx 

Ri 

if  E  =  0  Ri_n. 

For  n  =  1  (Ashida  and  Egashira  1977), 

dhb  =  3  +  \  dRi 

dx  Ri  3  dx 


For  mild  slopes,  uniform  flow  is  achieved  downstream,  Ri  becomes  for 
the  normal  depth  R  ;  and  dRi/dx  =  0  .  For  these  conditions,  Akiyama 
and  Stefan  (1981)  found 

~5  |3  +  f  +  ["(78,  0  +  +  4&S2tan  0  1 

R  =  ^  1  .  t _ LJ.J - 1 - £ - X  (44) 

n  2S2tan  G 

Since  R  is  independent  of  x  ,  Equation  (43)  can  be  integrated  with 
n 

respect  to  x  to  obtain: 

hb  =  f-  (x  '  Xo)  +  hbo=  EU  -  Xo>  +  hbo  (45) 

n 


where 

xq  =  reference  point  for  beginning  of  underflow 

h^  *  initial  depth  of  underflow  (In  practice,  the  initial  thickness 

of  the  underflow,  h,  ,  can  be  estimated  from  the  plunge 

bo 

depth,  hp  ,  see  Appendix  C.) 

From  Equation  (43),  h^  increases  linearly  downstream.  Equation  (44) 
shows  that  as  the  slope  increases  (i.e.,  R^  decreases)  and/or  entrain¬ 
ment  E  increases,  the  depth  of  the  underflow  increases  downstream.  If 
there  is  no  entrainment  (i.e.,  if  8=0),  then  the  depth  of  the  underflow 
remains  constant  at  h^.  The  underflow  velocity  and  density  at  point  x 
can  be  determined  from  Equations  (37),  (38)  and  (45). 
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103.  Hebbert  et  al.  (1979)  (see  also  Fischer  et  al.  1979  or 
Imberger  and  Patterson  1981)  derived  a  similar  relationship  for  a  tri¬ 
angular  cross  section  (Figure  21).  They  found 

hb  =  !  E(X  -  Xo)  +  bbo  (46) 


where  E  is  given  by  Equation  (35)  and 


5  tan  0  -  |  e  2p  u? 
_ j  _  a  b 

4E  +  5C,  /sin  <p  Ap,  gh, 
b  b  b 


(47) 


Equations  (35)  and  (47)  must  be  solved  iteratively,  when  C  ,  0  , 

3  b 

and  <t>  are  known  and  n  =  3.2  is  assumed.  The  dilution  AQ/q^ 

can  then  be  calculated  from 


M  = 

Qi 


(48) 


where  =  u^h^  tan  4>  for  triangular  cross  sections.  From  conserva¬ 
tion  of  mass 


(p  -  P.)Q.  =  constant 
a  i  i 


and  the  density  of  the  underflow  is  therefore 


P 

a 


(ca  ~  P1)Q1 

QA  +  aq 


(49) 


(50) 


104.  In  summary,  if  the  bed  slope  is  mild  (e.g.,  S  <  1/150), 
the  thickness  of  the  underflow  will  be  the  normal  depth  of  flow. 
Equations  (45),  (37),  and  (38)  can  be  used  to  predict  the  thickness 
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and  propagation  speed  of  an  underflow  in  a  rectangular  channel.  If  a 
triangular  cross  section  is  more  representative,  Equations  (48)  and 
(47)  should  be  used.  The  friction  factor  can  be  obtained  from  a  Moody 
diagram. 


Interflows  and  Intrusions 


105.  A  density  interflow  or  intrusion  occurs  when  a  density  cur¬ 
rent  leaves  the  river  bottom  and  propagates  horizontally  into  a  quies¬ 
cent  stratified  fluid.  The  density  current  enters  the  stratified  fluid 
at  the  level  where  the  densities  of  the  two  fluids  are  equal.  This 
fact  distinguishes  intrusions  from  overflows  and  underflows  in  which 
the  densities  of  current  and  ambient  fluid  are  not  the  same.  Once  the 
intrusion  enters  the  stratified  water  column,  all  turbulence  (e.g. , 
bottom  generated)  quickly  collapses  (i.e.,  dissipates)  and  the  in¬ 
trusion  assumes  the  properties  of  the  water  column.  The  density  of 

the  underflow  at  the  depth  where  the  density  current  enters  the  water 
column  and  becomes  an  intrusion  may  not  be  the  same  as  the  density  of 
the  river  water  entering  the  reservoir,  because  of  the  initial  mixing 
and  entrainment.  Changes  in  flow  rate  and  density  of  the  current  due 
to  initial  mixing  and  entrainment  must  therefore  be  known  before  an 
intrusion  can  be  analyzed. 

106.  Imberger  et  al.  (1976)  assumed  that,  since  the  intrusion  is 
neutrally  buoyant,  exact  symmetry  exists  between  the  inflow  and  outflow 
problem,  provided  the  inflow  momentum  is  small.  Using  the  dimensionless 
parameter  R  which  is  defined  as: 

R  =  F  Gr1/3  (51) 

m 

where 

2 

F  *  internal  Froude  number  =  q  /N  L 
m  2  A  2  rn  m 

Gr  =  Grashof  number  =  N  L  /« 

m  v 

q  =  inflow  per  unit  width  at  intrusion  takeoff  point  after  entrain- 
m  _ 

ment 
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cT  / 

S - -£•  =  \l& —  — S  =  buoyancy  frequency  averaged 

p  dz  W  p  h  .  '  .  , 

m  "mm  sion  thickness  h 

m 


over  intru- 


Ap  *  change  in  density  over  h 


m 


m 


density  at  centerline  of  intrusion 
average  vertical  eddy  viscosity 


and  the  Prandtl  number 


Pr 


\/KT 


Imberger  et  al.  (1976) 


where  k^,  *  average  vertical  diffusivity  of  heat, 
defined  three  flow  regimes: 

a.  R  >  1  .  The  intrusion  is  governed  by  a  balance  between  in¬ 
ertial  and  buoyancy  forces.  The  intrusion  length  L  is 

m 

proportional  to  time  t  such  that 


m 


0.44  (qN  )1/2t 


(52) 


The  intrusion  speed  L  /t  and  thickness  h  are  therefore 

m  m 


u  -  L  /t  =  0.44  (qN  )1/2  =  0 
mm  th  m 


1/2 


(53) 


and  substituting  in  u^  -  Q^/wt^  (assuming  uniform  flow)  and 
solving  for  h^  yields 


h  = 
m 


(54) 


b. 


—5 /6 

Pr  <  R  <  1  .  The  flow  regime  is  dominated  by  a  viscous- 
buoyancy  balance.  The  intrusion  length  is 


L 

m 


0.57LR 


2 / 3 t 5/6 


(55) 
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and  the  intrusion  thickness  is 


h 

m 


=  5.5L  Gr 
m 


-1/6 


(5;) 


In  this  regime,  64  percent  of  the  layer  lies  above  the  center- 
line  (Imberger  1980).  Therefore,  the  thickness  of  the  intru¬ 
sion  above  the  centerline  is 


h  !  =  3.5L  Gr~1/6 

m^u  m 


(57) 


and  below  the  centerline  is 


=  2.0L  Gr~1/6 
m 


(58) 


c. 


R  <  .  This  flow  regime  is  dominated  by  vicosity  and 

diffusion.  The  intrusion  length  is 


L  =  CLR3/4Gr“1/8 
m 


(59) 


and  the  thickness  is  of  the  order  LGr  3^R3^3.  The  coefficient 
C  is  unknown. 

107.  The  coefficients  in  Equations  (52)  and  (55)  were  obtained 
using  laboratory  data  from  Zuluaga-Angel  et  al.  (1972),  Maxworthy  (1972), 
and  Manins  (1976)  and  verified  with  data  from  Silvester  (1979).  The 
predictive  equations  were  also  compared  with  field  data  from  Cherokee, 
Fontana,  and  Wellington  Reservoirs.  Cherokee  and  Fontana  Reservoirs 
were  in  the  inertial  regime  (R  >  1),  while  Wellington  was  in  the  viscous 
regime.  According  to  Imberger  (1980),  Equations  (54),  (57)  and  (58) 
underpredict  the  thickness  of  the  flow  zone. 

108.  Manins  (1976)  obtained  similar  results  for  the  inertial 
range  and  went  on  to  show  that  internal  waves  (shear  waves) 
affect  the  intrusion  thickness  and  to  a  lesser  extent  the  intrusion 


66 


speed.  The  effect  of  the  internal  waves  is  to  adjust  the  isopycnals 
(lines  of  constant  density)  so  that  they  return  to  a  horizontal  position 
at  steady  state  (Mahony  and  Pritchard  1977). 

109.  Unless  dissolved  solids  dominate  the  density  profile  (i.e., 
high  Pr) ,  intrusions  in  most  reservoirs  fall  into  the  inertial  range  and 
Equations  (52),  (53),  and  (54)  should  be  used.  These  equations  were 
therefore  used  to  analyze  the  interflows  summarized  in  Table  1.  The 
intrusion  centerline  was  assumed  to  be  the  depth  in  the  lake  where  the 
lake  temperature  (i.e.,  density)  and  inflow  temperature  (density), 
including  the  effects  of  inflow  mixing,  matched.  The  total  inflow  per 
unit  width  was  obtained  by  dividing  the  total  inflow  (including  inflow 
mixing)  by  the  average  reservoir  width  of  the  zone  of  conveyance  at  the 
intrusion  centerline  level.  The  upper  and  lower  half-thicknesses  of  the 
intrusion  were  then  computed  separately  using  Equation  (34)  and  assuming 
the  total  inflow  per  unit  width  was  equally  divided  above  and  below  the 
centerline  (i.e.,  q^  was  divided  by  2  in  Equation  54  and  the  resulting 

.  h  was  the  half-thickness).  This  procedure  was  followed  to  account  for 
m 

temperature  (density)  profiles  that  were  not  symmetrical  about  the 

intrusion  centerline.  The  actual  application  of  Equation  (54)  involved 

assuming  a  half-thickness  h  ,  ,  calculating  the  Ap  over  h  ,  from 

nrs  m  m^ 

the  temperature  (density)  profile,  and  finally  calculating  h^  .  This 
process  was  continued  until  the  assumed  h^  converged  on  the  computed 
hm*s  *  which  normally  occurred  within  3  iterations  with  judicious  choices 


of 


110.  Initial  computations  using  the  maximum  flow  rate,  as  recom¬ 
mended  for  the  plunge  point  calculation,  resulted  in  intrusion  thick¬ 
nesses  that  were  too  large.  Subsequent  analyses  indicated  that  correct 
intrusion  thicknesses  could  be  computed  if  the  flow  rates  were  averaged 
over  the  time  scale  required  for  the  intrusion  to  pass  through  the  lake. 

This  time  scale  was  computed  from  Equation  (53)  with  L  equal  to  the 

m 

length  of  the  lake,  and  Ap  and  h  taken  across  the  metalimnion. 

m  m 

For  typical  summer  stratification  conditions  in  DeGray  Lake,  the  ap¬ 
propriate  time  scale  was  4  days.  During  the  spring  and  fall,  when  the 
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stratif ication  is  weaker,  the  time  scale  increases  to  6  or  more  days. 
When  a  time  scale  of  4  days  for  the  summer  months  and  6  days  for  early 
spring  and  late  fall  was  used,  the  following  results  were  obtained 
from  Equation  (54): 


Observed  _  Predicted 


Date 

h  .  ,  m 

m^u 

h  ;  ,  ,  m 
mil 

h  j  ,  m 
m*iu 

*W’  m 

31  Aug 

76 

1.8 

2.2 

1.8 

2.0 

17  Jun 

77 

3.0 

2.0 

2.8 

2.4 

14  Nov 

78 

3.8 

2.3 

3.4 

2.8 

1  May 

79 

5.5 

4.0 

5.8 

5.5 

16  Jul 

79 

2.4 

1.1 

2.4 

1.4 

23  Jul 

81 

2.0 

2.0 

1.9 

1.9 

An  example 

of 

these  computations  is  given 

in  Appendix  B. 

111. 

The 

total  intrusion 

thicknesses  predicted  by 

Equation  (54) 

were  usually  within  0.5  m  or  20%  of  the  observed  values,  provided  the 
correct  time  scale  was  used  to  average  the  storm  flows  and  the  width  of 
the  zone  of  conveyance  was  used  to  obtain  the  discharge  per  unit  width. 
For  any  specific  reservoir,  the  time  scale  can  be  computed  from  Equation 
(53)  assuming  equals  the  reservoir  length,  Ap^  equals  the  average 

density  difference  between  the  epilimnion  and  hypolimnion,  and  h^ 
equals  the  average  thickness  of  the  metalimnion.  Selection  of  the  width 
for  the  zone  of  conveyance  requires  extensive  knowledge  of  the  reservoir 
hydrodynamics  and  flow  patterns.  In  general,  side  arms,  and  dead-end 
coves  should  not  be  included  in  the  conveyance  width  (Figure  23). 

Aerial  overflights  after  major  storm  events  (e.g.,  Figure  5)  may  provide 
insight  into  flow  patterns.  For  hand  computations  of  specific  events, 
conveyance  widths  should  be  obtained  from  cross-sectional  data  in  the 
region  of  the  reservoir  of  interest  (e.g.,  the  upper  end).  If  the 
interflow  enters  the  pool  at  an  elevation  where  uncleared  timber  may 
determine  the  conveyance  limits,  sonar  soundings  such  as  those  reproduced 
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Figure  23.  Zone  of  conveyance  for  DeGray  Lake,  AR. 


in  Figure  6  may  also  be  helpful.  For  mathematical  model  simulations 
(e.g.,  CE-QUAL-R1) ,  the  conveyance  width  should  be  selected  to  be  repre¬ 
sentative  of  the  lower  region  of  the  reservoir  near  the  dam. 

112.  Once  the  intrusion  thicknesses  are  computed  from  Equation 
(54),  the  intrusion  speed  can  be  computed  from  Equation  (53)  using  the 

total  density  difference  Ap  across  the  total  intrusion  thickness 

m 

h  .  The  following  results  were  obtained: 
m 


Date 

31  Aug  76 
17  Jun  77 
1  May  79 
16  Jul  79 
23  Jul  81 


Observed 
Speed,  m/sec 
0.05 
0.06 
0.06 
0.03 
0.03 


Predicted 
Speed,  m/sec 
0.05 
0.06 
0.04 
0.05 
0.03 
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The  observed  speeds  were  for  movement  of  the  center  of  mass  and  not  the 
leading  edge,  which  moved  at  a  faster  rate.  For  the  storm  flows,  the 
observed  speed  was  averaged  over  3  or  4  days  and  did  not  consider  the 
period  when  the  flows  stalled  (paragraph  39) . 

113.  In  general,  the  predicted  intrusion  speeds  were  of  the  cor¬ 
rect  order  of  the  magnitude.  It  must  be  emphasized,  however,  that  the 
accuracy  of  these  predictions  is  highly  dependent  on  the  accuracy  of  the 
thickness  computations.  Equation  (53)  assumes  that  buoyancy  forces 
drive  the  intrusion.  Therefore,  in  reservoirs  dominated  by  advective 
forces,  the  accuracy  of  the  predictions  may  be  questionable. 

114.  Equations  (54)  and  (53)  are  recommended  tc  compute  intrusion 
thicknesses  and  speeds.  They  can  be  used  for  both  base  and  storm  flows 
provided  the  flows  are  averaged  over  the  proper  time  scale  and  the  width 
for  the  zone  of  conveyance  is  used. 
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PARI  IV:  SUMMARY  AND  RECOMMENDATIONS 


Review  Summary 

115.  A  literature  review  of  field  studies  on  density  currents  in 
reservoirs  produced  the  following  findings: 

a_.  Analysis  of  field  data  on  density  currents  is  complicated  by 
irregular  reservoir  morphometry,  unsteady  and  multiple  trib¬ 
utary  inflows,  variable  stratification,  and  man-made  pertur¬ 
bations. 

t).  Old  roadways,  railroad  crossings,  cofferdams,  causeways, 

constrictions,  and  other  hidden  obstructions  can  cause  irregu¬ 
larities  in  or  blocking  of  the  anticipated  flow  field, 
c^.  In  many  reservoirs,  the  flow  path  follows  the  old  thalweg, 
with  incomplete  lateral  mixing.  The  appropriate  width  and 
length  scales  therefore  appear  to  be  the  width  and  length  of 
the  old  river  channel  or  zone  of  conveyance. 
d.  Transport  into  and  out  of  coves,  embayments,  etc.,  is  enhanced 
with  pool-level  fluctuations.  When  the  pool  level  is  rising, 
water  moves  back  into  the  embayments.  During  periods  of  fal¬ 
ling  pool  levels,  water  moves  out  of  the  coves, 
je.  When  inflows  enter  a  wide  lake  (i.e. ,  width  >  1  km)  in  the 

northern  hemisphere,  Coriolis  accelerations  cause  the  density 
interfaces  to  tilt  to  the  right  side  of  the  basin. 
f_.  Inflowing  water  tends  to  pool  at  the  plunge  point,  which  can 
result  in  water  quality  problems.  The  pooling  enhances  delta 
formation. 

£.  The  location  of  the  plunge  point  is  highly  dynamic.  It  can 

move  several  kilometers  in  a  few  hours  in  response  to  dynamic 

flow  events  such  as  a  storm  event  or  hydropower  generation. 

-3 

h.  For  reservoirs  characterized  by  mild  slopes  (e.g.  10  ), 

initial  mixing  and  entrainment  are  small,  averaging  about  25 
percent . 
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jL.  Interflows  (intrusions)  stall  and  collapse  when  there  is  no 
hydraulic  gradient  (i.e.,  inflow  or  outflow). 
j_.  Even  when  a  stable  underflow  or  interflow  occurs,  inflow 

materials  can  still  be  transported  to  the  surface  waters  by 
wind-generated  mixing  mechanisms  and  changes  in  meteorological 
conditions . 


Predictive  Techniques 

116.  All  of  the  predictive  techniques  analyzed  in  this  study 
worked  well  under  laboratory  conditions  (i.e.,  steady  and  uniform  flow, 
uniform  channels  and  geometry,  and  linear  stratification);  however,  the 
application  of  these  techniques  to  existing  reservoirs  is  complicated 
by  irregular  morphometry,  varying  inflow  characteristics,  varying 
stratification,  and  reservoir  operation.  Many  reservoirs  are  dendritic 
in  shape  and  characterized  by  nonuniform  widths,  varying  bottom  rough¬ 
ness,  islands,  and  multiple  tributaries.  Inflow  characteristics,  such 
as  flow,  temperature,  and  TDS  and  SS  concentrations,  vary  diurnally, 
synoptically  (storm  events),  and  seasonally.  In-lake  density  strati¬ 
fication  is  dynamic,  nonlinear,  and  sometimes  unstable.  Changes  in 
reservoir  operation  (i.e.,  rule  curve,  withdrawal  depth,  release  rate, 
etc.)  can  have  a  significant  effect  on  inflow  density  currents  and  are 
not  considered  in  any  of  the  predictive  techniques. 

117.  The  recommended  predictive  techniques  are  described  hereafter 
and  sample  applications  are  presented  in  Appendix  B. 

Overflows 

118.  Overflows  occur  when  the  inflowing  water  density  is  less  than 
the  reservoir  water  surface  density.  This  condition  typically  occurs  in 
the  spring  when  rivers  warm  more  rapidly  than  reservoirs.  Predictive 
equations  are  available  for  the  separation  point  (Equation  1),  overflow 
thickness  (Equation  10),  and  frontal  propagation  speed  (Equation  8). 

The  separation  point  is  analogous  to  the  plunge  point.  Since  overflows 
are  exposed  to  wind-induced  mixing  mechanisms  and  heat  transfer  at  the 
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air-water  interface,  predictions  are  susceptible  to  large  errors  and 
should  be  used  with  caution.  Predictions  of  separation  point  depths  for 
the  overflows  in  Table  1  were  in  error  by  over  100  percent. 

Plunge  Point  Location 

119.  Density  inflows  plunge  beneath  the  water  surface  when  the 
buoyancy  forces  balance  the  advective  forces.  The  plunge  point  location 
is  determined  by  computing  the  hydraulic  depth  at  which  the  flow  plunges. 
Several  investigators  have  proposed  formulas  for  this  depth  asstiming 
steady  state,  uniform  velocity  and  density  profiles,  and  constant  reser¬ 
voir  density,  among  other  assumptions.  The  formulas  typically  depend  on 
the  inflow  densimetric  Froude  number,  bed  slope,  and  friction  coef¬ 
ficients.  Akiyama  and  Stefan  (1981)  (Equation  28)  included  initial 
mixing,  which  increased  the  depth  of  the  plunge  point.  Equations  (18), 
(22),  (24),  and  (28)  predicted  the  plunge  point  depths  to  within  +  1  m 
(i.e.,  1  km  along  the  longitudinal  axis)  for  both  base  flow  and  storm 
events.  No  specific  formulation  is  recommended  since  no  one  formula  is 
superior  for  all  reservoirs.  The  formula  of  Hebbert  et  al.  (1979) 
(Equation  22)  is  used  in  the  computer  algorithm  (Appendix  C) ,  since 

it  employs  a  triangular  cross  section  and  there  is  a  unique  relation¬ 
ship  between  reservoir  depth,  width  and  longitudinal  distance.  It  is 
recommended,  however,  that  the  width  of  the  zone  of  conveyance  be  used 
in  Equations  (18),  (24),  and  (28)  if  known. 

Underflows 

120.  After  the  inflow  plunges,  it  moves  along  the  submerged  river 
channel  as  an  underflow.  Turbulence  generated  at  the  bed  entrains 
ambient  reservoir  water  into  the  underflow.  If  the  bed  slope  is  mild 
(e.g.,  S  <  1/150),  Equations  (45),  (37),  and  (38)  can  be  used  to  predict 
the  thickness  and  propagation  speed  of  an  underflow  in  a  rectangular 
channel.  If  a  triangular  cross  section  is  used.  Equations  (46)  and  (48) 
can  be  used. 

Interflows  (Intrusions) 

121.  Intrusions  occur  when  an  underflow  leaves  the  old  river 
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channel  and  propagates  horizontally  Into  a  stratified  water  column  at  a 
level  of  neutral  buoyancy.  The  intrusion  problem  is  therefore  analogous 
to  the  withdrawal  problem  in  a  stratified  fluid.  Equations  (54)  and 
(53)  can  be  used  to  estimate  the  intrusion  thickness  and  speed,  respec¬ 
tively,  for  the  nonviscous  regime.  The  width  for  the  zone  of  conveyance 
should  be  used  to  compute  the  discharge  per  unit  width,  and  storm  flows 
should  be  averaged  over  the  appropriate  time  scale,  which  is  computed 
from  Equation  (52). 
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APPENDIX  A:  DENSITY  OF  WATER 


1.  The  density  of  water  is  a  function  of  temperature,  total  dis¬ 
solved  solids,  and  suspended  solids: 

pw  "  PT  +  ApTDS  +  APSS  (A1) 


where 


Ap 

Ap 


TDS 

SS 


3 

water  density,  kg/m 

3 

water  density  resulting  from  temperature  variations,  kg/m 

3 

density  increment  due  to  total  dissolved  solids,  kg/m 

3 

density  increment  due  to  suspended  solids,  kg/m 


2.  Many  different  formulations  have  been  proposed  for  determining 
the  density  variation  due  to  temperature.  Several  mathematical  models 
used 


PT  =  1000. 


(T  -  3. 98) 2  (T  +  283) 
(503.57)  (T  +  67.26) 


(A2) 


where 

O- 

T  *  water  temperature,  C. 

Eq.  (A2)  attains  a  maximum  density  of  1000.  kg/m3  at  T  *  3.98°C. 

3.  A  simpler  and  computationally  more  efficient  relationship  is 

PT  -  999.841  +  6.59583xl0"2T  -  8.45123xl0'3T2  (A3) 

+  5.29159  x  10"5T3 

Chen  and  Millero  (1977)  proposed 

PT  -  999.8395  +  6.7914x10"2T  -  9.0894xl0-3T2  (A4) 

+  1.0171xl0_4T3  -  1.2846x10"6T4  +  1. 1592xlO-8T5 


A1 


-  5.0125x10  L1T6 


4.  The  density  increment  due  to  total  dissolved  solids  is 

ApTDS  =  PTDS  *  CTDS 

where 

3  3 

PTDS  =  density  conversion  factor,  kg/m  /g/m 

3 

CTDS  =  TDS  concentration>  g/m 

Reported  values  for  PTDg  vary  from  less  than  0.0005  to  0.001;  a  typical 
value  is  0.00078.  PTDS  will  depend  on  the  chemical  makeup  of  the  TDS 
and  temperature.  Chen  and  Miller  found 

PTD$  =  8.221xl0-4  -  3.87xl0"6T  +  4.99xlO~8T2  (A6) 

5.  The  density  increment  due  to  suspended  solids  is 

Apss  *  Css  (1  -  1/SG)  x  10‘3  (A7) 

where 

3 

CSS  =  88  concentration,  g/m 

SG  =  specific  gravity  of  the  suspended  solids 
For  a  SG  =  2.65,  Eq.  (A7)  reduces  to 

Apss  =  0.00062  Css  (A8) 

6.  Chen  and  Millero  also  modified  the  density  for  pressure  (i.e., 
at  depth) : 

P___  -  P  /(I  -  P/K)  (A9) 

wp  w 


A2 


APPENDIX  B:  EXAMPLE  COMPUTATIONS 


1.  In  this  appendix,  two  examples  from  DeGray  Lake,  Arkansas,  are 
presented  to  illustrate  the  use  of  the  recommended  predictive  equations 
for  plunge  point  location,  initial  mixing,  interflow  thickness,  and  pro¬ 
pagation  speed.  The  first  example  is  a  dye  study  which  was  performed  in 

3 

July  1979  under  steady  flow  conditions  (Q^  =  4.6  m  /sec),  where,  after 

plunging,  the  dye  moved  into  the  reservoir  as  an  interflow  (Figure  Bl) . 

The  second  example  is  an  unsteady  storm  event  which  occurred  in  June 

3  3 

1977,  when  the  flow  increased  from  approximately  5  m  /sec  to  370  m  /sec 
in  8  hours  (Figure  B2)  and  the  turbid  storm  water  entered  the  lake  as  an 
interflow.  Both  of  these  examples  are  fully  described  in  Johnson 
et.  al.  (1983). 


July  1979  Dye  Study 

2.  At  the  time  of  dye  injection,  the  following  conditions  were 
measured : 

3 

Flow  rate  =  4.6  m  /sec 

Inflow  temperature  T  =  27.9°C 

1  2 

River  cross-sectional  area  A  =  7.2  m 

River  top  width  W  =  24.7  m 

Hydraulic  depth  =  A/W  =  0.3  m 

Lake  surface  temperature  T^  =  29.5°C 

Using  Equation  (A2)  for  the  density  of  water, 

”  996.288  kg/m3 

p  =  995.733  kg/m3 
a 

The  inflow  densimetric  Froude  number  (Equation  2)  is: 

F  ,  _Jli _ (A. 6/7. 2) _ 

1  (^g  hi)Jj  J-(5  6  x  10“4)*9. 8*o.  3]^ 

where 

ui  "  VA 


Bl 


Figure  Bl.  Dye  concentrations  in  ppb  for  the  July  1979 
dye  study,  DeGray  Lake,  Arkansas. 
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Figure  B2.  June  1977  storm  event,  DeGray  Lake,  Arkansas 
(Percent  T  is  percent  transmittance). 
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Api  996.288  -  995.733  ,  ,  ,„-4 

ei  =  p - 9957733 -  5,6  x  10 

a 


3.  The  plunge  point  depth  is  computed  from  Equation  (18)  (Savage 
and  Brimberg  1975) : 


h 


P 


1/3 


where 


F 

P 


2.05 

(l+«) 


ar 


2.05  (  0.0012  \  0,478 

(1+0.2)  \  0.02  / 


0.45 


Therefore, 


h 


P 


(4.6/70)' 


l 9. 8*5. 6  x  10~4*(0.45)2 


1.6  m 


The  width  of  the  zone  of  conveyance  was  70  m, which  was  the  bank-full 
width  of  the  old  river  channel  (Figure  6) .  The  observed  hydraulic 
plunge  point  depth  was  1.8  m. 

4.  The  initial  mixing  coefficient  is  computed  from  Equation  (32) 
and  the  downstream  densimetric  Froude  number  from  Equation  (44) 


2S2S 

isS1P  +  ft  +  +ft)2+  48S2s]  ** 


where  -  0.25,  S2  *  0.75,  and  g  «  0.0015  are  site-independent 
coefficients.  Therefore,  for  any  specific  reservoir  the  mixing  coef¬ 
ficient  and  downstream  densimetric  Froude  number  are  constant  such  that 


B4 


5.  The  centerline  depth  of  the  interflow  is  determined  by  com¬ 
paring  the  inflow  temperature  (density)  with  a  measured  temperature 
(density)  profile  in  the  lake  (Figure  B3) .  The  centerline  depths  for 

Y  -  0  ,  y  -  0.03  ,  and  y  =  0.10  are  5.9,  5.7,  and  5.5  m,  respective¬ 
ly.  The  observed  centerline  depth  was  approximately  5.5  m.  It  should 
be  noted,  however,  that  the  inflow  temperature  varied  from  27.4°  to 
28.4°C  during  the  period  of  dye  injection.  This  range  corresponds  to 
the  depth  range  5.0  to  6.1  m  which  falls  within  the  inflow  zone 
(Figure  B3) . 

6.  Using  a  centerline  depth  of  5.5  n  (y  ■  0.10),  the  upper  and 
lower  limits  of  the  inflow  are  determined  from  Equation  (54)  assuming 
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DEPTH 


Figure  B3.  Comparison  of  predicted  and  observed  interflow  placement  and 
thickness,  July  1979  dye  study,  DeGray  Lake,  Arkansas. 
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one-half  of  the  flow  is  distributed  above  the  centerline  and  one-half 
below  the  centerline  depth.  Therefore, 


h 


m%u 


(Qbo/(2W))2  \1/3 
g  Ap  ] 

T*  / 


where  Q.  =  5.06  m  /sec  and  W  =  70  m  (i.e.,  the  width  of  the  zone  of 
bo 

conveyance).  Since  Ap^  is  dependent  on  h^^,  this  equation  is 
solved  by  first  selecting  an  h^^,  then  determining  Ap^  from  Figure 


B3  and  solving  for  h 


m^u 


This  procedure  is  continued  until  the  cal- 

For  example,  initially 


culated  h  .  converges  on  the  initial  h  . 

m^u  ro>su 

select  h  .  =  2.1  m  .  The  temperature  at  z  *  5.5  -  h  ,  =3. 4m  is 

m^u  -  dfsu 

29.1°C,  and  p(29.1°)  =  995.941  kg/m  .  Therefore,  Apy/p  =  (996.202  - 

995. 941)/996. 202  =  2.62  x  10_4  and  h,  =  2.4  m.  Since  2.4  m  is  not 

nrsu 

similar  to  2.1  m,  the  computation  is  repeated  using  h  ,  =  2.3  m  .  The 

temperature  and  density  at  the  upper  limit  are  now  29.2°C  and  995.912, 

-4 

respectively.  Therefore  Ap  /p  =  2.91  x  10  and  h  ,  =  2.3  m.  The 

u  m^u 

lower  limit  is  similarly  computed.  First,  select  h^  ^  =  0.9  m,  which 
corresponds  to  a  temperature  of  26.2°C  and  density  of  996.756  kg/m^. 

Since  the  computed  h^^  is  1.9  m  and  not  similar  to  the  initial  h^^ 
of  0.9  m,  the  procedure  is  repeated  with  h^^  =  1-.4  m  .  This  half¬ 
thickness  corresponds  to  a  temperature  of  24.0°C  and  density  of  997.323 
3 

kg/m  .  Solving  for  h^^  results  in  1.5  m,  which  is  similar  to  the  1.4  m 
initially  selected.  The  total  thickness  of  the  interflow  is  1.4  m  + 

2.4  m  *  3.8  m  ,  which  agrees  with  the  vertical  dye  distribution  (Figure 
B3) . 

7.  The  speed  at  which  the  interflow  propagates  into  the  reservoir 
is  given  by  Equation  (53): 


u 

m 


0.194 
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where  Ap  and  h  are  over  the  total  thickness  of  the  interflow. 
m  3  3 

Therefore,  Ap  =  1.41  kg/m  ,  p  =  996.202  kg/m  ,  h  =  3.8  m  and  u  = 

m  mm 

O.O45  m/s  =4.5  cm/s  .  The  leading  edge  of  the  dye  plume  was  observed 
to  advance  at  3.3  cm/s. 


June  1977  Storm  Event 


8.  All  of  the  predictive  equations  reviewed  in  this  report  assumed 
steady  state.  In  this  example,  they  are  applied  to  an  unsteady  storm 
event  (Figure  B2)  to  determine  their  accuracy.  The  required  data  are: 


Flow  rate 
Inflow  temperature 


Q  =  370  m  /a 


19. 9°C 


Inflow  density  p.(19.9°)  =  998.251  kg/m 

T  =  28.2°C 
a 


Ambient  lake  temperature 

a 

Ambient  lake  density  pa<28.2)  *  996.202  kg/mJ 
Other  coefficients  and  parameters  which  are  site-specific  and/or  in¬ 
variant  include: 


Total  friction  coefficient  f  =  0.024 
Bottom  slope  S  *  1.2  x  10 
Conveyance  width  W  *  440  m 
Profile  constant  =  0.25 
Profile  constant  ■  0.75 
Entrainment  coefficient  8  =  0.0015 
9.  In  this  example,  the  plunge  depth  is  computed  using  the  equa¬ 
tions  of  Savage  and  Brimberg  (1975)  (Equation  18)  and  Hebbert  et  al. 
(1979)  (Equation  20).  Since  the  densimetric  Froude  number  F^  used  in 
the  Savage  and  Brimberg  equation  is  invariant  for  a  specific  reservoir, 
Fp  *  0.45  from  paragraph  3.  Therefore,  the  Savage  and  Brimberg  plunge 
depth  is 


/ I2  V'3 . 

/  (370/440)2  \  1/3 

Uiv 

\(9. 8*0. 002*0. 452)/ 

5.5  m 


where 


Api  998.251  -  996.202 

ci  “  p  996.202 

a 


0.002 


The  Hebbert  et  al.  formulation  assumes  a  triangular  cross  section  and 
the  site-specific  densimetric  Froude  number  is  given  by  Equation  (19) : 

u 

„2  sin  S  tan  (1  -  0.85C,/2  sin  S) 

*  -  r  u 

n  CD 


where 

CQ  =  total  drag  coefficient. 

Substituting  the  site-specific  characteristics  for  DeGray  Lake  yields 


j~sin(0. 0012) 

0.055 


tan (85°) 


(1-0. 85(0. 055)'2sin(0. 0012) 


.5 


and  the  equation  for  the  plunge  depth  (Equation  20)  becomes 


3 

Substituting  =  0.002  and  =  370  m  /sec  yields  13.3  m,  which  must 

be  divided  by  2  to  obtain  the  hydraulic  depth  for  a  triangular  cross 
section.  Therefore,  h^  =  13.3/2  *  6.7  m  ,  which  is  comparable  to  the 
5.5  m  computed  from  Savage  and  Brimberg  and  the  observed  hydraulic  depth 
of  5.4  m.  These  two  computations  indicate  the  differences  in  magnitude 
that  can  be  obtained  from  the  different  formulations  for  the  plunge 
depth.  It  is  also  worth  noting  that  the  densimetric  Froude  number 
required  for  both  formulations  differ  slightly  (i.e.,  0.45  versus  0.50). 

10.  Since  the  downstream  densimetric  Froude  number  computed  from 
Equation  (44)  is  site  specific  and  not  directly  a  function  of  flow,  the 
computed  initial  mixing  coefficient  for  this  example  will  be  the  same  as 
the  previous  example  (i.e.,  y  ■  0.03).  This  value  compares  with  a 
measured  mixing  coefficient  of  0.07. 
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11.  Assuming  no  initial  mixing  and  an  inflow  temperature  of 
19.9°C,  the  centerline  of  the  intrusion  would  be  at  an  elevation  of 
116.7  m  (Figure  B4) .  With  an  initial  mixing  coefficient  of  0.07  and  a 
surface  temperature  of  29°C,  the  revised  inflow  temperature  would  be 
20.5°C  and  the  centerline  of  the  intrusion  would  be  at  117.0  m.  As 
shown  in  Figure  B2,  the  inflow  temperatures  varied  from  approximately 
19°C  to  24°C  during  the  period  of  elevated  flow  and,  assuming  no  mixing, 
would  have  entered  the  reservoir  as  separate  intrusions  between  the 
elevations  of  116.3  and  119.1  m.  This  zone  closely  approximates  the 
intrusion  thickness  (Figure  B4).  For  the  intrusion  computations,  a 
centerline  at  117.0  m  will  be  assumed. 

12.  The  flow  rate  required  for  the  intrusion  thickness  computation 
(i.e..  Equation  54)  must  be  averaged  over  the  time  scale  given  by  Equa¬ 
tion  (53) : 


t 

s 


«s 


where 

L  *  reservoir  length 

Ap  *  density  change  across  metalimnion 
m 

h  *  average  thickness  of  metalimnion 
m 

For  DeGray  Lake,  under  average  summer  stratification  conditions, 
L  *  32,000  m  ,  *  8  m  ,  Ap^  *  p(25°C)  -  3  kg/m^  ,  and 


t 

s 


32,000 


3.4  x  10  sec 


0.194 


1000 


*9.8*8 


3.9  days 


3 

The  resultant  average  flow  rate  for  period  16-19  June  was  63  m  /sec. 

13.  The  upper  and  lower  half-thicknesses  are  determined  from 
Equation  (54),  using  a  centerline  elevation  of  117  m  and  an  average 

3 

flow  rate  of  63  m  /sec.  For  the  upper  half-thicknesses  h  ,  first 
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assume  h^u  *  3  m.  From  the  temperature  profile  on  19  June  (Figure 

B4),  Ap  /p  =  (p(20.5°)  -  p(26.2°))/p(20.5°)  =  (998.126  -  996.756)/ 
m  ~ 

998.126  =  1.37  x  10  .  Therefore, 


m^u 


/  (63/440/2)2  \  1/3  _  2 
3,0  V9. 8*0. 00137  /  "  2,2 


Since  2.2  <  3.0  m 

,o 


profile  at  25.3  C 
1.14  x  10-3,  and  h 

i 

less  than  the  assumed 
In  this  case. 


try  =  2.6  m,  which  intersects  the  temperature 

Therefore,  Ap  /p  *  (p(20.5°)  -  p(25.3u))/p(20.5°)  = 
tn 


m^u 


=  2.3  m  .  Since  the  predicted  h 


ra^u 


m^u 


(  i .  e .  >  2.3  <  2.6  m)  ,  try  h 


Ap  /p  =  (p (20. 5°)  -  p(24.7°))/p(20.5°)  =  9  x‘l0-4 
m 


is  still 
=  2.4  m. 
and 


*\n*4u  converges  to  2.5  m. 


For  the  lower  half-thicknesses,  first  assume 
Then,  from  Figure  B4, 


h  ,  ,  «■  2.5  m. 
m*sl  o  -3 

p(20.5°)  ■  1.06  x  10  ,  and  h^^  =  2.4  m, which  is  similar  to  the 


Ap  /p  =  (p(14.6°)  -  p(20.5°) ) / 
m 


assumed  h^  of  2.5  m.  The  computed  intrusion  zone  is  shown  in  Figure 
B4. 

14.  The  intrusion  speed  is  computed  from  Equation  (53) 


VAPm 

—  gh 
p  m 


using  h^ 
4.9  m  and 


*\nSsu  +  ^m^sl 


and  Ap  across 


Ap  /p 

m 


1.96  x  10 


-3 


m 

yields 


m 


h^  .  Substituting  h^  * 

*  0.06  m/sec,  which  is  the 


same  as  the  observed  speed  of  0.06  m/sec. 


Discussion 


15.  The  numbers  used  in  the  two  preceding  examples  differ  slightly 
from  those  used  in  the  main  text  of  this  report  because  different  in¬ 
lake  conditions  were  selected.  Figure  B4  clearly  illustrates  the  vari¬ 
ations  that  can  be  expected  in  temperature  and  intrusion  characteristics. 
It  is  therefore  important  to  use  in  the  calculations  inflow  and  in-lake 
characteristics  that  were  collected  at  the  same  time  and  location. 


2 


Because  of  the  dynamic  nature  of  these  characteristics,  differences  in 
predictions  can  be  expected. 


APPENDIX  C:  COMPUTER  ALGORITHM 


1.  This  appendix  outlines  a  mathematical  inflow  algorithm  to  be 
used  in  a  one-dimensional  water  quality  model  such  as  CE-QUAL-R1  (Envir¬ 
onmental  Laboratory  1982)  to  predict  inflow  placement,  plunge  point 
depth,  initial  mixing,  underflow  entrainment,  and  intrusion  thickness. 
The  algorithm  is  based  on  the  formulations  reconmended  in  the  main  text 
of  this  report  and  considers  a  reservoir  to  be  composed  of  n  hori¬ 
zontal  layers  of  variable  thickness  Az(k)  and  variable  density  p(k) 
where  k  =*  l,n  .  The  layers  are  numbered  from  the  bottom  (k  •  1)  to  the 
water  surface  (k  *  n) .  The  algorithm  requires  the  following  information 
for  each  tributary: 

3 

a.  Inflow  rate  Q  ,  m  /sec . 

1  3 

t>.  Inflow  density  ,  kg/m  .  (This  should  be  a  function  of 
temperature,  TDS,  and  SS.  See  Appendix  A.) 

£.  Conveyance  width  W(k)  for  each  horizontal  layer  k  ,  m 
(should  be  representative  of  main  body  of  reservoir) . 

d.  Bottom  slope  S  ,  m/m  (computed  from  reservoir  depth  Z  ,  and 

—  m 

length  L  along  the  old  thalweg:  S  *  Z  /L  ). 

tn 

£.  Half  angle  <J>  representative  of  the  upper  end  of  the  reser¬ 
voir  (see  Figure  21) . 

f.  Total  friction  coefficient  for  underflow  f  (f  *  f.  +  f.), 

—  t  t  b  i 

Upper  mixed  layer  depth  zmix>  m  (composed  of  layers, 
k  =  m,n  ) . 

2.  For  each  tributary,  the  following  computational  steps  are 
required: 


a.  Compute  the  plunge  point  depth  h^  using  the  formulation  of 
Hebbert  et  al.  (1979),  (Equation  20): 


where 


(Cl) 


Cl 


„  2  sin  S  tan  4>  d  -  0. 85CL*  sin  S) 
F  .  “  U 

p  CD 


and  £i  ‘  "“[(“a  * 


p  «  reservoir  water  surface  density,  kg/m 

CD  ’  V4 

Determine  the  layers  in  the  plunge  zone  (i.e.,  k  *  p,  n  where 
p  is  the  layer  at  depth  h  ).. 

Ik  Compute  the  initial  mixing  coefficient,  y  ,  using  the  formula¬ 
tion  of  Jirka  and  Watanabe  (1980)  (Equation  32)  and  the  den- 
simetric  Froude  number  computed  from  Equation  (C2): 


Y  =*  1.2  F  -  0.2 
P 


£.  Compute  the  revised  rate  and  density  pbo  after  initial 

mixing: 


Q±  (1  +  Y) 


Tl>.  +  -i 

<1  +  Y> 


If  Pbo  <  pfl  ,  a  density  overflow  occurs. 

(1)  If  h  >  Z  ,  (i.e.,  p  <  m),  distribute  the  inflow 

p  mix 

throughout  the  depth  Z  ^  (k  -  m,n) ,  assuming  a 
maximum  velocity  at  the  water  surface  and  a  lower  inflow 
zone  limit  m  .  No  initial  mixing  occurs  since  it  is 
assuned  Fp  is  small  (<  1)  and  the  density  gradient  at  the 
bottom  of  the  mixed  layer  is  sufficiently  strong  to 
inhibit  mixing  and  force  the  inflow  into  the  mixed  layer 
only. 


C2 


(2)  If  hp  <  (i.e.  p  >  m)  ,  distribute  the  inflow  rate 

after  initial  mixing  QbQ  throughout  the  depth  h 
(k  *  p,n) ,  assuming  a  maximum  velocity  at  the  water 
surface  and  a  lower  inflow  zone  limit  p  (Figure  Cl). 


Figure  C2.  Definition  diagram  for  an  underflow. 


underflow  reaches  a  point  of  neutral  buoyancy  (i.e.,  = 

p(k))  or  until  the  dam  is  reached  (k  *  1)  and  the  inflow 
is  inserted  in  the  bottom  layer.  For  any  layer  k  (see 
Figure  C2)  the  depth  at  the  end  of  the  layer  is 

Vk>  '  +  hbo(k)  <C9) 


Since  volume  and  mass  are  conserved, 


Vk> 


(^doVg 


(pbQ(k)  -  p(k)) 
p(k) 


and 


Pb(k> 


pbo(k)Qbo(k)  +  P(k>Qe(k) 

Qb(k) 


where  the  entrained  flow  is 


(CIO) 


(Cll) 


C5 


Qe(k)  -  Qb(k)  -  Qbo(k) 


(C12) 


Equations  (C9) ,  (CIO),  (Cll) ,  and  (C12)  are  recomputed 
for  each  layer  until  (k)  <  p(k)  ,  then  an  intrusion 
occurs  with  layer  k  as  the  layer  of  neutral  buoyancy. 
_f.  The  thickness  of  intrusion  is  determined  from  Equation  (54). 
Since  the  vertical  density  gradient  in  the  vicinity  of  an 
intrusion  is  usually  nonlinear  (Figure  C3) ,  the  intrusion 


Figure  C3.  Definition  diagram  for  an  interflow. 


thickness  computation  is  separated  for  the  upper  and  lower 
limits.  The  upper  half-thickness  is  therefore: 


*Vsu 


3.0 


(Qbo(k)/(2W(k)))2 

p(k) 


1/3 


(C13) 


C6 


where 


k  =  number  of  layer  that  intrusion  enters 

Q,  (k)  =  total  inflow  rate  with  all  entrainment,  but 
bo  ’  ^ 

averaged  over  the  appropriate  time  scale,  m  / sec 
W(k)  =  width  of  layer  k,  m 

Ap  =  density  difference  over  the  half-thickness 
u  3 

h  Ull  ,  kg/m 

Q,  (k)  must  be  divided  by  2  in  Equation  (C13)  since  it  is 
bo 

assumed  that  half  the  flow  is  above  the  level  of  neutral 
buoyancy  and  half  below.  Q^Ck)  must  also  be  averaged  over 
the  time  scale 


0.194^^jn  ghm  ^  ^ 


where  Ap  is  the  change  in  density  across  a  metalimnion  of 
m 

thickness  h 

m 


Equation  (C13)  is  solved  with 


Apu  =  p(k+i)  -  p(k) 


m 

h^u  =  ^  (Az  (k-l+i)  +  Az(k+i))/2 
i=l 


adding  one  layer  at  a  time  until 


ABS 


3.0 


(k)/(2W(k)))2 

rAp 


Az 


where  Az  =  average  layer  thickness,  m 

The  lower  limit  is  determined  in  a  similar  manner. 

The  interflows  and  overflows  are  distributed  between  upper 
and  lower  limits  according  to  parabolic  velocity  profiles 
used  in  selective  withdrawal. 
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APPENDIX  D:  NOTATION 


A 

B 

C. 


SS 

“TDS 


cross-sectional  area,  m 

2 

top  width  of  the  underflow,  m 
bottom  drag  coefficient 
total  drag  coefficient 

coefficient  to  account  for  nonhydrostatic  pressure  distribution 

3 

SS  concentration,  g/m 

3 

TDS  concentration,  g/m 
entrainment  coefficient 

uk 

underflow  densimetric  Froude  number 


(£b*V 


=  inflow  densimetric  Froude  number  = 


=  internal  densimetric  Froude  number 


t 

Gr 

g 


VNmL 


normal  depth  densimetric  Froude  number 
densimetric  Froude  number  at  plunge  point 
bed  friction  coefficient 
interfacial  friction  coefficient 
total  friction  factor  =  f,  +  f . 

D  1 

2  4  2 

Grashof  number  *  N  L  /< 
m  v 

2  . 

gravitational  acceleration,  m  /sec 

Pb  '  pa 

relative  acceleration  of  gravity  *  - 


g,  m  /sec 


2 

"b 

bbo 


observed  plunge  depth,  m 
observed  depth  of  underflow,  m 
depth  (thickness)  of  underflow,  m 
initial  depth  of  underflow,  m 


critical  depth 


2  1/3 


m 


D1 


hi 

hm 

*\n*sl 

^m^u 

h 

n 

h 

P 

h 

s 

K 

L 

L 

m 

L 

s 

m 

N 

P 

N 

c 

N2 


P 

Pr 

Q 

Qb 

V 

<>e 

<>i. 

QHOE(k) 

qi 

Sn 

R 


inflow  depth,  m 
intrusion  thickness,  m 
lower  intrusion  thickness,  m 
upper  intrusion  thickness,  m 
downstream  normal  depth,  m 
plunge  point  (separation)  depth,  m 

thickness  of  the  underflowing/overflowing  current,  m 

constant  *  ft/S2 

length  of  reservoir,  m 

intrusion  length,  m 

overflow  length,  m 

lower  uptake  zone  limit 

dimensionless  ratio  *  h  /h 

P  n 

dimensionless  ratio  =  h  /h 

c  n 

buoyancy  frequency  averaged  over  intrusion  thickness 

v  -2 

h  ,  sec 
m 

pressure,  bars  (1  bar  -  10.2  m  of  water) 

Prandtl  number  *  kv/kt 

3 

discharge,  m  /sec 

3 

underflow  discharge,  m  /sec 

revised  flow  rate  including  entrainment,  downstream  of  the 

3 

plunge  point,  m  / sec 

3 

rate  of  entrainment  flow,  m  /sec 

3 

original  inflow  rate,  m  /sec 

3 

entrained  inflow  from  layer  k,  m  /sec 

3 

unit  width  discharge,  m  /m/sec 

inflow  per  unit  width  at  takeoff  point  after  entrainment, 
m/m/sec 

1/3 

dimensionless  parameter  ■  F  Gr 

n 

normal  depth  Richardson  number 


D2 


g,hb 

Ri  =  Richardson  number  =  — 

% 

r  =  entrance  mixing  ratio 

S  =  bottom  slope 

SG  =  specific  gravity  of  the  suspended  solids 
Sl*  ^2  =  profile  constants 

o 

T  =  water  temperature,  C 

=  lake  surface  temperature,  °C 
=  inflow  temperature,  °C 
=  inflow  temperature  after  mixing,  °C 
t  =  time,  sec 

u^  =  underflow  velocity,  m /sec 

u^o  =  initial  underflow  velocity,  m/sec 

u^  =  inflow  velocity  =  q^/h^,  m/sec 

u  =  intrusion  speed  =  L  /t,  m/sec 

m  m 

Up  =  velocity  at  the  plunge  point,  m/sec 

u  =  propagation  speed  of  the  overflow  front,  m/sec 

s 

W  =  width,  m 

W  =  width  of  zone  of  conveyance,  m 

c 

W  =  initial  width,  m 

o 

W  =  bank-full  width  of  old  river  channel,  m 

r 

W  =  width  at  distance  x,  m 

x 

x  *  distance  from  source,  m 

x  =  reference  point  for  beginning  of  underflow,  m 

o 

z  =  reservoir  depth,  m 

Z  =  maximum  reservoir  depth,  m 

m 

^mix  *  uPPer  mixed  layer,  m 

a  =  ratio  of  interfacial  to  bed  friction  *  f,/f 

i  o 

g  *  coefficient  in  entrainment  equation 

Y  *  initial  mixing  coefficient 

Az  *  average  layer  thickness,  m 

Az(k)  -  thickness  of  layer  k,  m 


D3 


f 


Ap 


m 


Ap 

Ap 


SS 

TDS 


e 


b 


e 


i 


p(k) 


p(z) 

pa 

pb 

pbo 


density  difference  of  ambient  and  underflow,  kg/m 

density  difference  between  receiving  water  and 

3 

inflow  «  (p  -  p  ),  kg/m 

3  X 

3 

change  in  density  over  h  ,  kg/m 

m  3 

density  increment  due  to  suspended  solids,  kg/m 

3 

density  increment  due  to  total  dissolved  solids,  kg/m 
relative  density  difference 


relative  density  difference 


Ap,  to 

b  a 

Ap  ,/p 

i  a 


average  vertical  diffusivity  of  heat,  m  /sec 

2 

average  vertical  eddy  viscosity,  m  /sec 
ratio  of  normal  depth  to  plunge  point  depth  ■ 
mixing  efficiency 
bed  angle 


h  /h 
n  p 


a/(l  +  a) 

3 

density  of  layer  k,  kg/m 

3 

density  variation  with  depth,  kg/m 

3 

ambient  density  of  the  reservoir  surface  waters,  kg/m 

3 

density  of  the  underflow,  kg/m 

revised  inflow  density  after  entrainment  at  the  plunge 

3 

point,  kg/m 

3 

entrainment  density,  kg/m 

3 

density  of  inflowing  water,  kg/m 

3 

density  of  intrusion,  kg/m 

3 

density  of  reservoir  surface  waters,  kg/m 

■} 

water  density  resulting  from  temperature  variations,  kg/m 

33 

density  conversion  factor,  kg/m  /g/m 

3 

water  density,  kg/m 

half-angle  of  the  triangular  cross  section 


V 


D4 


